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i.o Introduction 


1.1 Welcome 

This is the User Guide to the CIRS instrument of the Cassini spacecraft. A User Guide is 
being prepared for each of the Cassini instruments by the instrument science teams. We 
hope you find this guide useful, and that is enables you to conduct research and publish 
findings based on the Cassini CIRS public dataset. Enjoy! 

1.2 Purpose 

This guide is written for scientists who are using, or contemplating the use of publicly 
available CIRS data in the Atmosphere’s Node of the Planetary Data System (PDS). It is 
intended to pass along the ‘corporate knowledge’ of CIRS-team insiders to the wider 
scientific community, especially in the longer term when the Cassini mission has ended and 
the instrument teams disbanded. Specific goals of this document include: 

• An introduction to the instrument, and the types of observations made for 
different target bodies and various science goals. 

• What data is available, and how to extract the particular subset you are 
interested in. 

• Examples of how the data can be used. 

• Pitfalls and traps, especially calibration uncertainties and noise interferences in 
the data. 

This guide for the CIRS instrument is one of several guides available, for various instruments 
on the Cassini orbiter and Huygens probe. 


1.3 Scope 

This guide focuses on CIRS data: its acquisition; organization, including ancillary 
information such as pointing, and metadata; some details of data processing 
(cleaning/calibration); extraction of data, especially binary data; and some examples of 
usage. However this document is not intended to be a detailed technical guide to the 
instrument itself, nor a mathematical or computer algorithm guide to the details of the 
Fourier Transform and calibration process. 


1.4 How to use this guide 

This user guide may be either read end-to-end as a complete introduction to the CIRS 
PDS dataset, or used as a reference guide for the experienced user to consult regarding 
details of specific topics. 
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2.0 General Information 


This chapter first describes the Cassini Mission in general and the CIRS instrument in 
particular. That is followed by a detailed account of the various scientific objectives and 
measurements made by CIRS, and then we tackle the topics of data gaps and instrument 
performance, including interferences and other artifacts that affect the data. 

2.1 Instrument Characteristics 

This section describes the Cassini mission and CIRS instrument. 

2.1.1 Cassini Mission 

The Composite Infrared Spectrometer is a dual-interferometer carried on the Cassini 
spacecraft Remote Sensing Palette (RSP). Cassini was launched on 15th October 1997 and 
made four gravity assist maneuvers before reaching Saturn: Venus (twice), Earth and finally 
Jupiter. The slow, distant Jupiter flyby (137 RJ) at the end of 2000 allowed for a six-month 
observing campaign that was used to evaluate the performance of the Cassini instruments, 
and achieve some early science. The CIRS Jovian dataset forms the first part of the CIRS 
archive. 

After achieving Saturn orbit on July 1 st 2004, Cassini commenced a four-year ‘Prime 
Mission’ lasting until July 2008. It included 76 spacecraft orbits (or revolutions: ‘revs’) of 
Saturn known as 0, A, B, C, 3 ... 76. The original orbits 1 and 2 were replaced with A, B and C 
when the mission was redesigned in order to accommodate radio link issues with the 
Huygens probe, discovered after launch. The probe delivery to Titan occurred on rev C. 
There were 44 targeted (i.e. close) fly-bys of Titan during the first 78 orbits. 


Table i: Cassini Mission Phases 


Phase 

Dates 

Orbits 

(revs) 

Sequences 

Titan 

Flybys 

Other 

Encounters 

Prime 

07/01/04-06/30/08 

A-78 

S1-S41 

T0-T44 

See Figure 

Extended (EM) 

07/01/08-06/30/10 

78-133 

S42-S60 

T45-T70 

See Figure 

Solstice (SM) 

06/30/10- 

134-270 

S61-S99 

T71-T126 

See Figure 


The second mission phase was the Equinox Mission (EM), sometimes called the 
Extended Mission (XM), over 2 years from July 2008 to July 2010. This mission phase was 
designed to reach the Saturn northern spring equinox in August 2009, and comprised a 
further 55 orbits and 26 further flybys of Titan (T45 to T70). At the time of writing, the third, 
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and final anticipated mission phase is the Solstice Mission (SM, CSM, or XXM for second 
extended mission), beginning in July 2010 and continuing through to the anticipated 
spacecraft disposal in Saturn’s atmosphere in mid-2017. This phase will encompass the 
northern summer solstice of Saturn in early 2017, and includes 136 orbits, 56 further Titan 
encounters, and multiple encounters with other moons (see Table 1 and Figure 2 below). The 
last part of the mission features ‘proximal orbits’, where the spacecraft periapse (closest 
approach to Saturn) will be moved from outside to inside the main rings, providing a very 
close look at Saturn’s atmosphere. 




Cassini Mission Overview 





Four-Year Prime Tour, Equinox Mission, and Solstice Mission, May 2004 - September 2017 
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Figure 2: Cassini mission overview showing Saturn orbits and satellite encounters, as of April 22 nd 2012. 


TWTs and OSTs 

The time allocation between various Cassini instruments (also known as the time-line 
‘integration’ process, as opposed to the later ‘implementation’ process, which developed 
actual pointing designs) was carried out by a number of virtual teams, known as either TWTs 
(Target Working Teams) or OSTs (Orbiter Science Teams). These were: 

• TOST - the Titan Orbiter Science Team 
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• SOST - the Satellites Orbiter Science Team 

• XD TWT - the Cross-Discipline Target Working Team 

• Saturn TWT - the Saturn Target Working Team 

• Rings TWT - the Rings Target Working Team 

• MAG TWT - the Magnetospheres Target Working Team 

CIRS maintained representation on all of these virtual (not physically co-located) working 
teams that met by teleconference at regular intervals. Following timeline integration, it was 
the responsibility of each team to develop spacecraft pointing commands for the allocated 
time-period, rather than by request to a central design office (as on the Galileo Mission) a 
system known as ‘Distributed Operations’ (Distributed Ops). 

For further information on the Cassini mission, see the MISSION. CAT file in any of the 
Cassini PDS volumes. 

2.1.2 Instrument Design 

CIRS consists of two interferometers; a far-IR polarizing interferometer sensitive from 10- 
600 cm 1 (1000-17 pm), and a mid-IR interferometer sensitive from 600-1400 cm ' (17-7 pm) 
(see Figure 3 and Figure 4 below). These two share common components including the 
telescope (Cassegrain, with 508 mm primary beryllium mirror, and 76 mm secondary), fore- 
optics, reference laser and scan mechanism, which moves the retroreflectors of both 
interferometers simultaneously. Engineering details are given in (Kunde et al., 1996). 
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Figure 3: CIRS mechanical layout from Kunde et al. (1996). 

The far-IR radiation is sensed by a large (3.9 mrad FWHM) thermopile detector (actually 
two, one for reflected and one for transmitted beam at the polarizer/analyzer), known as 
FPi. The mid-IR radiation falls on two 1X10 pixel FlgCdTe arrays, known as focal planes 3 & 4. 
FP3 is an array of photoconductive (PC) detectors sensitive from 600-1100 cm 1 (17-9 pm), 
and FP4 is an array of photovoltaic (PV) detectors sensitive from 1100-1400 cm 1 (9-7 pm). 
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There is no ‘FP2’ as it was dropped from the original design. Further details of the mid-IR 
arrays can be found in (Nixon, Infrared Limb Sounding of Titan with Cassini CIRS, 2009). 



Figure 4: schematic layout of the CIRS instrument optical pathways. Incoming light from the telescope (not shown) at top 
is split into the mid-IR (left path) and far-IR (right path) interferometers. Both interferometers share a common scan 
mechanism. 


Each of the 1X10 mid-IR arrays (Figure 5) has a 5-channel amplifier/signal processor, 
meaning that only 5 outputs from the 10 detectors in each array can be recorded at a time. 
The maximum number of simultaneous interferograms from CIRS is therefore 11 (1 from FPi, 
5 from FP4, 5 from FP4). Flowever, there are observing modes in which some of the focal 
planes are not used, leading to a lesser number of active channels. 
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Figure 5: CIRS mid-infrared focal plane 
layout, showing detector spacing and 
numbering. 


Sampling is triggered by fringes of a reference red 
laser (784.510.3) nm and the raw interferogram is later 
decimated by a numerical filter (factors 18, 4.5, 5.0 for 
FPi, 3, 4 respectively). A white light source is used to 
sense the zero-path-difference (ZPD) point when both 
arms of the interferometer are equal length. Scanning 
occurs over a maximum path difference of about 2 cm, 
but may be less. Spectral resolution is determined by the 
mirror travel distance. Velocity is constant, 1 so this 
equates to scan time, in units of 1/8 second (1 RTI, or 
‘real time interrupt’, from the spacecraft clock). 2 The 
resolution is varied from approximately 15.5 cm-i (lowest 
resolution, for a -5 sec scan, or 40 RTI) to 0.5 cm-i 
(highest resolution, from a 50 s scan, 400 RTI). 3 

See the INST.CAT file on the archive volumes for 
additional details. 


2.1.3 Instrument operation modes 


Set-point and thermal stability 

While most of the instrument operates at -170 K, stabilized by heaters, in contrast the 
mid-IR focal plane is cooled to 75-85 K by a passive radiative cooler, which is linked to the 
thermally-isolated mid-infrared focal plane assembly (FPA) by a ‘cold finger’. This operating 
temperature is controlled by replacement heaters commanded to different set-points (see 
Table 2) beginning with SP 22 (76.22 K) and rising in increments of 0.4 K. Note that the actual 
temperatures achieved by the instrument at the various set-points differ from the 
commanded value: for example at set-point 22 CIRS actually stabilizes at 76.301. 


1 Nominally at least. Practically, it is feedback-controlled to within +/- 1%. 

2 Really, the scan time is controlled in Flight Software to the millisecond level. As of FSW v6, 
5 ms is added to each scan to help control spike artifacts, see “2.5 Electrical interferences 
(spikes) and other data issues”. 

3 Maximum allowed range is 16-416 RTI (FSW VI-V5) and 24-416 RTI (FSW v6 on). 
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Table 2: CIRS set point temperatures (all in Kelvin) 


Set 

Point # 

FPA 

Temp 

Min 

Temp 

Max 

Temp 

Set 

Point # 

FPA 

Temp 

Min 

Temp 

Max 

Temp 

31 

69.52 

69.32 

69.72 

15 

81.37 

81.17 

81.57 

30 

70.27 

70.07 

70.47 

14 

82.10 

81.90 

82.30 

29 

71.02 

70.82 

71.22 

13 

82.83 

82.63 

83.03 

28 

71-77 

71-57 

71-97 

12 

83.56 

83.36 

83.76 

27 

72.52 

72.32 

72.72 

11 

84.28 

84.08 

84.48 

26 

73-26 

73.06 

73-46 

10 

85.01 

84.81 

85.21 

25 

74.00 

73.80 

74.20 

09 

85-74 

85-54 

85.94 

24 

74-75 

74-55 

74-95 

08 

86.47 

86.27 

86.67 

23 

75-49 

75-29 

75-69 

07 

87.19 

86.99 

87-39 

22 

76.22 

76.02 

76.42 

06 

87.92 

87.72 

88.12 

21 

76.96 

76.76 

77-16 

05 

88.64 

88.44 

88.84 

20 

77-70 

77 - 5 ° 

77-90 

04 

89-37 

89.17 

89-57 

19 

78.43 

78.23 

78.63 

03 

90.10 

89.90 

90.30 

18 

79-17 

78.97 

79-37 

02 

90.90 

90.70 

91.10 

17 

79.90 

79-70 

80.10 

01 

91.60 

91.40 

91.80 

16 

80.63 

80.43 

80.83 

00 

92.30 

92.10 

92.50 


The set-point is usually set to be as low as possible for best detector performance, 
subject to the expected heating load on the passive radiative cooler. For example, during a 
‘rolling downlink’ (spacecraft pointed with -Z to Earth, high gain antenna (HCA) to Earth, 
but rolling about the Z-axis for fields and particles science) the cooler may experience small, 
cyclic heating flux from a warm body (Sun, Saturn, or rings) in excess of its capacity to shed. 
Therefore, during these events, CIRS may either not take data, or take data but with an 
elevated set-point, commanded to be above the maximum temperature excursion predicted 
(see figure). 
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[INSERT: figure of thermal predict during rolling downlink, before and after set-point 
management] 

Shutter 

The instrument may command the shutter closed (for calibration) or open (default). The 
shutter affects only the mid-infrared detectors (FP3 and FP4), not the far-infrared detector 
(FPi). The shutter temperature is close to the stabilized instrument temperature (-170 K). 
Shutter open and closed modes appear to cause different baselines (‘shape functions’) in 
the FPi interferograms for electrical reasons (or possibly straylight?), although the optical 
path has not been occluded. See “5.1.5 Field of view modeling” for further details. 


FP3 and FP4 Pixel Modes 

Due to the fact that fewer amplification channels are available (5) than the numbers of 
detectors on each focal plane (10), a choice must be made about which detectors are in use 
an any given time. The choices are ‘ODD’, ‘EVEN’, ‘CENTERS’ and ‘PAIRS’. In the ‘ODD’ mode, 
odd numbered detectors are used (1,3, 5, 7, 9) and in the even mode, even-numbered 
detectors are used (2,4,6,8,10). As the numbering of the FP3 and FP4 detectors are reversed, 
FP3 detector 1 is side-by-side with FP4 detector 10. Therefore, to match spatial positioning, 
FP3 EVEN is usually selected along with FP4 ODD, and vice versa. In ‘CENTERS’ mode 
detectors 3-7 are used from FP3 and detectors 4-8 are used from FP4, which are again co- 
located. In ‘PAIRS’ mode, the signals from each pair of detectors within each array - 1812, 
3&4, 5&6, 7&8, 9&10 - are co-added before being sent into the amplifier, so that the arrays 
function as five-element arrays with 2x larger rectangular pixels but increased S/N sensitivity. 


Table 3: pixel modes 


Virtual 

Detector 

Number 

Focal 

Plane 

fp 3 

Physical 

Detector 

FP3 Mode 

Virtual 

Detector 

Number 

Focal 

Plane 

FP4 

Physical 

Detector 

FP4 

Mode 

00 

1 

1 






01 

3 

18(2 

PAIR 

21 

4 

1 

ODD 

02 

3 

3&4 

PAIR 

22 

4 

2 

EVEN 

03 

3 

5&6 

PAIR 

23 

4 

3 

ODD 

04 

3 

7&8 

PAIR 

24 

4 

4 

EVEN 

05 

3 

98(10 

PAIR 

25 

4 

5 

ODD 
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06 

3 

3 

CENTER 

26 

4 

6 

EVEN 

07 

3 

4 

CENTER 

27 

4 

7 

ODD 

08 

3 

5 

CENTER 

28 

4 

8 

EVEN 

09 

3 

6 

CENTER 

29 

4 

9 

ODD 

10 

3 

7 

CENTER 

30 

4 

10 

EVEN 

11 

3 

1 

ODD 

31 

4 

4 

CENTER 

12 

3 

2 

EVEN 

32 

4 

5 

CENTER 

13 

3 

3 

ODD 

33 

4 

6 

CENTER 

14 

3 

4 

EVEN 

34 

4 

7 

CENTER 

15 

3 

5 

ODD 

35 

4 

8 

CENTER 

16 

3 

6 

EVEN 

36 

4 

1&2 

PAIR 

17 

3 

7 

ODD 

37 

4 

3&4 

PAIR 

18 

3 

8 

EVEN 

38 

4 

5&6 

PAIR 

19 

3 

9 

ODD 

39 

4 

78(8 

PAIR 

20 

3 

10 

EVEN 

40 

4 

9800 

PAIR 


Co-add mode 

The co-add mode is either ‘on’ or ‘off’. If co-add is ‘on’ then an interferogram is stored in 
a buffer, and co-added to the next interferogram from the same detector, before being 
passed along the data chain. This reduces effective temporal sampling, but reduces the data 
volume by a factor 2 when data volume limitations are imposed. The length of an 
interferogram from two coadded scans corresponds to the length of the shorter of the two 
scans being coadded. The data from a coadded interferogram includes the frame scan 
diagnostics data block and frame scan ID block information from both of the original scans, 
but only one frame scan data block containing the coadded data. 

2.1.4 Types of Data at PDS 

The data selected for submission to the PDS falls into four broad categories. 
Housekeeping tables (HSK) are included to provide insight into the health and status of the 
instrument. Interferogram-related tables (DIAG, FRV, IFGM, IHSK, OBS) provide the 
interferograms themselves plus a wide-range of parameters describing the instrument 
modes, conditions, and settings that occurred during the scan. The navigational tables (GEO, 
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POI, RIN, TAR) detail the positions of the major bodies, describe where the instrument was 
pointed, and what bodies were within the FOV. Finally, the ISPM tables provide the 
calibrated spectra resampled onto standardized wavenumber grids. This information is 
summarized in Table 4 below. 


Table 4: PDS directories 


PDS Volume Directory 

Description 

Data Table Types 

DATA/HSK_DATA 

Raw housekeeping data at 64 sec 
intervals. 

HSK 

DATA/UNCALIBR 

Unprocessed raw interferograms 
(IFGM), basic instrument parameters 
(IFISK, OBS, FRV) and preliminary 
assessment of quality (DIAC). 

IFGM, IFISK, OBS, DIAG, 
FRV 

DATA/NAV_DATA 

Navigation and pointing data derived 
from the Cassini SPICE kernels and 
computed for the Cl RS observations. 

GEO, POI, RIN, TAR 

DATA/APODSPEC 

Calibrated and apodized spectra 

ISPM 


Descriptions of the fields within each table as well as the binary format of the tables are 
given in the DATASIS.PDF document available at any of the CIRS PDS node sites. 

2.1.5 Data Unit Definition 

The following table lists the types of data units that can be found in the DATA directory. 
These are delivered to the PDS, and archive at the PDS Atmospheres Node, in variable-length 
binary form. Flowever the tables have been translated into a mixture of ASCII and fixed- 
length binary records for increased accessibility, and are available in this form at the PDS 
Rings Node. 

Table 5: CIRS PDS data object types 

FISK These tables store the data from the housekeeping packets (e.g., 

instrument temperatures, voltages, and currents) that are required for the 
calibration of IFMs. Note that the housekeeping packets are created every 
64 seconds, and records do not have a direct one-to-one or one-to-many 
correspondence to science scans, as all other table records do. A full 
description of this data is in the CIRS flight software documentation. 

DIAC These tables contain diagnostic information determined regarding the raw 

interferograms, including the presence or absence of noise interferences 
(e.g., spikes, velocity variations, RWA), the identification of the laser mode, 
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and the location of the ZPD. 


FRV 

IFGM 

IHSK 

OBS 

GEO 

POI 

RIN 

TAR 


These tables store the voltages of the central fringe of the CIRS white-light 
interferometer, sampled at i second intervals after mirror flyback. 

These tables store the actual IFM records. Note that there will typically be 
eleven I FM records (one per signal channel) for a given SCET time. 

These tables store interpolated housekeeping data for use in the 
calibration. A small subset of fields from the FISK housekeeping tables are 
interpolated to the times of the actual scans. 

These tables store information relating to the instrument data-taking 
settings for each scan (e.g., RTI, co-add, shutter status) and the values of 
warning and indicator flags. 

These tables store the geometry, the relative positions and orientations of 
the spacecraft, primary body, and object bodies (which may include the 
primary, e.g. Saturn), usually at the time of the zero path difference (ZPD) 
of the interferogram (when most of the signal is collected). One entry is 
generated for every body in the target list. 

These tables store the pointing for each active detector at the ZPD of the 
IFM and at scan's end. One POI entry is generated for each detector and 
each non-rings non-primary target in the field of view. An entry is always 
generated for the primary as the target, whether or not it is in the field of 
view. Rings are described in a separate RIN table rather than as entries in 
the POI tables. 

These tables contain the pointing information for the rings. One entry is 
generated for each detector at all times. For the purpose of this table, the 
rings are assumed to fill the equatorial plane of the primary (z = o in the 
coordinate system of the primary with z along the rotation axis) out to a 
selected cutoff radius. If a line of sight intersects the primary body before 
intersecting the rings, or if the line of sight falls outside the cutoff radius, 
the values for most fields for that line of sight are set to -200. 

These tables describe which, if any, of the various targets are in the field of 
view of a given detector. Individual fields (e.g. TITAN) are set to 1 or 0 
according to whether the body appears in the FOV or not. Any, all, or none 
of these fields may be specified individually to select data with a particular 
subset of bodies in the FOV. Additionally, the field FOV TARGETS contains 
the same information as the individual byte fields, except stored as bits in a 
single 4-byte integer. This may provide a useful shorthand for selecting an 
exact target body combination. 
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ISPM 


These tables contain the interpolated (resampled) calibrated spectra. The 
older SPM type created formerly used a natural wavelength index resulting 
from the instrument, which meant that the final spectrum fell on irrational 
wavenumber values such 577.3256, 577.7792, instead of 577.250, 577,500, 
577-750 etc. The ISPM type contains the same spectra, but interpolated and 
re-sampled onto a more user-friendly grid. There is no loss of resolution or 
data in this process. NOTE: the point spacing of the ISPM tables does not 
correspond to the spectral resolution, nor is it half the resolution (Nyquist 
spacing). In practice the over-sampling is always smaller than Vi the 
spectral resolution. The ‘real’ resolution is contained in a table field 
(Nyquist, Rayleigh or FWHM, depending on which definition you prefer). 
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2.2 Science Objectives 

For a detailed description of the science objectives, see (Flasar, 2004). A short 
description is included here for guidance. The Cassini/Fluygens mission to the Saturnian 
system is designed to investigate the following targets: 

1. Saturn : Sensing of tropospheric and stratospheric temperatures, dynamics, and 
composition. This includes abundances of the major and minor species, the hunt for 
new gaseous species, the ortho-to-para ratio for hydrogen, isotope ratios for major 
species. Allocating time for Saturn observations was the responsibility of the Saturn 
TWT (Target Working Team), which had orbit (rev) and segment allocations divided 
up through negotiation with other working groups. The Atmospheres Working Group 
(AWG) also provided input and dealt with high-level science recommendations and 
mission objectives. 

2. Rings : Measurement of the infrared spectrum at varying phases, leading to 
conclusions about ring particle size, shape, composition, distribution and dynamics. 
Allocation of time for rings observing was carried out by the Rings TWT, in 
conjunction with the Rings Working Group (RWG). 

3. Titan : Sensing of tropospheric and stratospheric temperatures and composition. This 
includes abundances of the major and minor species, the hunt for new gaseous 
species, isotope ratios for major species, and dynamics. CIRS js^also able to sense the 
surface near sso cm 1 . Allocation of time for Titan observations was primarily done in 
the TOST group (Titan Orbiter Science Team), in conjunction with recommendations 
from the AWG. 

4. Other ('icy') satellites : Surface mapping in the infrared provides information on the 
surface composition, temperature, morphology, and age. Close passes of icy 
satellites resulted in time being allocated to the SOST (Satellites Orbiter Science 
Team), which in turn divided the time between teams. This effort was done with the 
oversight and guidance of the Surfaces Working Group (SWG). 

The overall mission may also be broken in physical target types: 

1. Atmospheres (Saturn and Titan) 

2. Surfaces (Titan, and other satellites). 

3. Magnetosphere: fields, particles and solar wind interaction (Sun-Saturn) 

4. Orbital mechanics (rings, satellites especially Flyperion). 


CIRS contributes to the investigation of many of these areas. 


2.3 Observation Strategies 

Here follows a description of how the various science investigations are implemented in 
terms of actual spacecraft pointing. Accompanying figures are provided in Appendix A. 
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2.3.1 Saturn 

CIRS Saturn atmospheric science goals are met with observations at a variety of 
inclinations, ranges from Saturn and, in some cases, in conjunction with other teams' 
observations. The standard CIRS requests include: 


Table 6: CIRS Saturn observation types 


Name 

Range 

(Rs=i 

Saturn 

radius) 

Details 

COMPSIT 

All ranges 
of R s 

Composition sit and stare at 0.5 cm' 1 spectral resolution with 
FP3/4 paired detector mode using either nadir viewing for low 
emission angle studies, or FP3 or FPi offset to obtain high 
emission angles. Uses all 3 focal planes. Search for oxygen 
compounds (C 0 2 and Fl 2 0 ) and new molecules. 

MIRMAP 

25-50 R s 

Mid-IR maps, centered at a single latitude, with the focal planes 
oriented N-S on Saturn to allow the planet's rotation to map 
out a latitude band. Uses all three focal planes and 3 cm 1 
resolution using FP3/4 paired detector mode. Arrays are 
centered on the central meridian. 

MIRTMAP 

25-50 R s 

Mid-IR maps that step across multiple latitudes with the focal 
planes oriented N-S on Saturn. Hybrid of MIRMAP’s to allow 
mapping temperatures across the entire planet regularly. 
Latitudes from 75°N to 75°S are covered with approximately 
2hrs of longitude coverage per latitude pointing. Uses all three 
focal planes and 3 cm' 1 resolution using FP3/4 paired detectors. 
Arrays always centered on the central meridian at each latitude 
position. 

FIRMAP 

15-30 R s 

Far-IR maps, constructed by slewing equator-ward over a single 
hemisphere. Slews overlap and repeat to map out the entire 
hemisphere. 15 cm 1 spectral resolution and FP3/4 paired 
detector mode. Some short temporal spacing between 
FI RM APs of the same hemisphere to search for thermal waves. 

FT RACK 

< 10 R s 

Campaign with VIMS and ISS to track features across the disk 
from limb to limb, including limb sounding. FPi nadir sounding 
at 3 cm' 1 , FP3/4 limb sounding at 15 cm' 1 . FP3/4 paired and blink 
detector modes. The goal is to study vertical structure and 
pole-ward heat transport. Only performed in the prime 
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LIMBMAP 


LIMBINT 


NADIROCC 


OCCLIMB 


RECMAP / 
COMPFT 


TEMPSIT 


MIRCMPSIT 


FIRCMPSIT 


mission. 

< 6 Rs Vertical sounding at multiple latitudes on the limb at 15 cm 1 . 

FP3/4 blink detector mode. Focal planes are placed normal to 
the local latitude limb to sound multiple altitudes. All 3 focal 
planes are used. Sometimes are labeled LIMBINT, if repurposed 
during the design phase. 

< 6 Rs Long integrations for vertical composition at one location on 

the limb, with 2 or 3 altitude positions. A spectral resolution of 
0.5 cm 1 until mid-2008 and 1 cm 1 thereafter. FP3/4 blink 
detector mode. 

< 10 R s Helium abundance measurements taken by following the Radio 

Science Team's occultation points across the disk. 
Observations performed at 3 latitude positions displaced by -2, 
0, +2 mrad along the central meridian longitude (CML) to allow 
for pointing uncertainties. Using 3 cm 1 spectral resolution and 
FP3/4 paired detector mode. 

< 6 Rs Independent verification of vertical temperature profile at 

Radio Science's occultation latitudes. Observations performed 
at 0.5 and 2.5 mrad above the 1 bar level (100 and 1 mbar?). 15 
cm 1 spectral resolution and FP3/4 blink detector mode. 

< 10 Rs Regional mapping and/or composition feature tracks, to fill in 

coverage not obtained by the other observations (0.5 or 3 cm 1 
spectral resolution and blink detector mode). Primarily used in 
the prime mission. 

Saturn approach temperature map. 3 cm 1 spectral resolution 
and paired detector mode. Focal planes oriented N-S on the 
central meridian. Planet rotates to cover all longitudes. Arrays 
cover Southern Hemisphere (north blocked by rings). 

Saturn approach composition sit and stare. 0.5 cm- 1 spectral 
resolution and paired detector mode. Focal planes oriented N-S 
on the central meridian. Arrays cover Southern Hemisphere 
(north blocked by rings). 

Saturn approach far-IR composition. 0.5 cm 1 spectral resolution 
and paired detector mode. FPi centered on South Pole and on 
mid-Southern latitudes. 
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2.3.2 Rings 

The name conventions for the ring observations have changed during the mission, driven 
partially by the team discovering that certain objectives were beyond the sensitivity of the 
instrument, and partially by the unanticipated degree of variation in effective temperature 
and emissivity with illumination conditions and observation geometry. With the exception of 
a few shifts of policy in observation planning that were made during main mission, most of 
the new naming conventions commenced with XM, and a few more with XXM. 

The science objectives are discussed individually in the first sub-section, below, and the main 
events driving the name changes that took place during the mission are discussed in 
afterwards. Table 10 presents a summary for how the naming conventions associated with 
the science objectives have changed, and 

-Table 9 summarizes the science objectives by observation name. 

Science Objectives and Observing Strategies 

Faint Ring Long Integrations: the low optical depths of the faint rings were taken with 15 
cm 1 integrations using FPi. The FMOVIE and FMONITOR observations taken during the Main 
Mission detected the F Ring only when averaging over many hours. The original geometry 
planned for nearly edge-on viewing, however it was discovered that, due to the thickness of 
the F ring, viewing from an inclination of 1.5 to 2 degrees maximized signal by allowing the 
ring to fill substantial portions of the FPi FOV. No opportunities were found for F Ring 
observations in XM, with the XXM still pending as of the time of this report. 

Composition Integrations: These observations were designed to search for absorption 
features between 50 and 1000 microns. Three observation types were planned. Long 
integrations at 0.5 cm 1 would target each of the main rings at ansa, from 20-40 Rs using FPi 
and FP3; FP3 measurements would take long integrations of the rings against Saturn, so that 
absorption features might be catalogued; and certain radial scans, “submillimeter” scans, 
would be made at 1 cm 1 spectral resolution of the lit and unlit sides of the rings to map the 
thermal characteristics and composition of the ring particles out to 1 mm. 

Stellar Occultations: Provide a direct measurement of the IR opacity of the rings. The only 
two targets that have been identified as being observable by CIRS are CW Leo (FP3) and Eta 
Carinae (FP4). These were planned for 15 cm 1 spectral resolution, centering the star in a 
single pixel of the array. 

Ring Thermal Properties: Constrain the overall thermal budget of the rings and how it 
drives the emission seasonally and during orbital periods. The flux into the rings depends 
upon the Solar elevation, distance from Saturn, and azimuthal location at points along the 
orbit, and the emission of the rings should vary in step with this. This variation will be 
catalogued using radial scans typically executed between 5 and 20 Rs over a range of 
spacecraft inclinations, from the lowest (5 0 ) to the highest possible inclination (75 0 ). Radial 
mapping (FPi, FP3) of the rings, on both lit and unlit sides, over a range of spacecraft 
elevations, inclinations, local times and phase angles, is performed to obtain broadband 
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radiometric measurements of the total flux in the CIRS wavelength range. Sets of 
observations are obtained in each of the inclined orbit intervals to map the temperature 
variation in the rings with changing solar illumination. Two types of scans are planned. So 
called temperature scans will consist of spectra at 15 cm 1 spectral resolution of the lit and 
unlit sides of the rings at many incidence and emission angles and provide prime information 
on the ring thermal gradient as a function of radial distance to Saturn. 

Ring Particle Thermal and Dynamical Properties: Infer surface properties, vertical dynamics 
and spin of ring particles. Observations of the cooling and heating of the ring particles 
entering and emerging from the planetary shadow are planned to derive particle thermal 
inertias for all three main rings. It will make measurements at moderate radial resolution 
(typically 1000 km) across the shadow boundaries at low spectral resolution (15 cm 1 ) with 
the FPi field of view. To constrain the vertical dynamics of ring particles, the temperatures of 
the main rings will be measured by CIRS along the ring azimuth of the main rings, from the 
exit of the shadow (morning) to the evening ansa, both on the unlit and unlit faces. This 
unique experiment will be realized with spectra at low spectral resolution (15 cm 1 ). Spins 
create both an azimuthal asymmetry in the ring temperature and a dependence of the 
temperature with the emission angle, due to day/night contrast. Circumferential scans at a 
variety of phase and emission angles will be executed to detect azimuthal asymmetries and 
the anisotropy in the ring particle emission function which are both function of particles spin 
and thermal inertia. Occasionally, when observing time is highly disputed, long azimuthal 
scans (8-to-20 hours long depending on geometry and face) will be replaced by a series of 
radial scans at different azimuths. 

Observation Types 

A summary of the observation names associated with the Prime Mission science objectives is 
listed in 


Table 7: summary of CIRS Prime Mission rings observations and scientific objectives 


Name: 

Scan Type: 

Objectives 

COMP 

Stare at A, B, or C ring ansa i cm 1 
resolution. 

Detect bands due to 
contaminants 

SUBM 

Radial scans at ansae, i cm 1 resolution, 
slow scans. 

Scan slowly at high resolution 
to map radial contaminants. 

TEMP 

Radial scans iscm 1 resolution 

Map radial variations of 
temperature. 

SHAD 

Azimuthal scans at fixed radii (usually 
R=82000 km, 108000 km, 132000 km), 15 

Observe azimuthal heating 
and cooling to constrain 
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cm 1 resolution. 

thermal inertia. 

VERT 

Azimuthal scans at fixed radii, planned 
in pairs for lit and unlit comparison. 15 
cm 1 resolution. 

Observe azimuthal heating 
and cooling on both sides of 
thick rings to find evidence 
for vertical transport of 
particles. 

FMOVIE 

Fixed spot on F ring for duration of orbit 

Observe azimuthal variation 
and dumpiness. 

FMONITOR 

Shorter version of F Movie 


CWLEOOCC 

Fix on RA/DEC of target star as it occults 
rings 

Determine optical thickness 
of rings in IR 

ZEROPHAS 

Map several locations in B ring to 
assemble phase curves near 0 phase 

A one-off follow-up 
discovery mode which failed 
to detect . 


During the mission the team learned the following lessons that changed the observation 
strategies during design implementation as follows: 

(i) Long integrations showed no hint of bands from contaminants; some SUBM 
radial scans became TEMP scans; fewer COMP scans were scheduled (except 
when, due to geometry, only the ansae were available). 

(ii) Noise limitations at low wavenumbers were prohibitive for observing the change 
from emission to reflection spectrum in the 5-20 cm 1 range. This resulted in 
further curtailment of the SU BM scans. 

(iii) The effective temperature and emissivity have much more directional variation 
and detailed dependence on illumination conditions, radius, and resolution than 
had been anticipated; in contrast, the azimuthal curves could be over-fit, for 
cooling and thermal properties. Consequently, at the onset of XM it was decided 
that fine radial resolution scans contained more information than fine azimuthal 
resolution at a given radius. The TMAP observations were started as thermal 
mapping campaigns to take multiple radial scans so as to provide more 
information on radial structure and use less time for over-fit azimuthal curves. 

(iv) The optimal viewing geometry for faint rings is not edge-on, but at a slight 
elevation so that the ring fills the FOV. The spectra from the original FMOVIES 
had to be averaged over whole orbital periods in order to get a signal at all, so 
more of these were not scheduled. 
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(v) Some observations suggested a thermal phase surge. One observation, a 
ZEROPHAS, was created to catalogue this peak at the end of Prime Mission, but 
obtained a null result and was not repeated. 

(vi) The phase behavior of the thermal emission and filling factor came to be of 
interest for fitting comprehensive models of ring emission. Consequently, in the 
XXM, observations are labeled by the phase/latitude bin that they cover. The 
[N,S]PXXLXX observations were designed to complete phase curve studies 
started at the end of the PM. The PIE observations in this category were 
specifically designed to capture the phase variation at Solar elevations where the 
spacecraft was in equatorial orbit during MM (i.e., for which there was no rings 
data). 

(vii) The competition for time around Saturn Equinox prompted the design of the 
EQ[N,S][L,S]B observations, short for EQuinox North/South Long/Short baseline. 
These were designed to detect the overall cooling rates of the rings on timescales 
of hundreds of days (Long Baseline) around equinox, and hours or days (Short 
Baseline) around it. They are radial scans at 15 cm 1 scheduled at geometries as 
self-similar as tour constraints allow. 

(viii) TDIFF measurements are designed in pairs to constrain vertical transport by 
observing radial scans at the same Hour Angle from the lit- and unlit sides of the 
rings at similar spacecraft elevations, solar elevations and phase angle. While 
vertical transport was an objective of the azimuthal VERT scans, the unanticipated 
variation in thermal emission with illumination and observation geometry 
resulted in these being somewhat ineffectual because the different geometries 
often resulted in a multiplicity of potentially influential parameters varying 
between the lit- and unlit observations. 

(ix) SHAD[UN,LIT] were introduced in XXM to make up for a lack of high-latitude 
radial scans of the shadow region during PM and XM. 

Summaries of the observation names associated with the science objectives for XM and 
XXM are in Error! Reference source not found, and Table 9. A tabulation of how the science 
objectives and observation names changed during the full mission is in Table 10. 


Table 8: summary of CIRS XM Mission rings observations and objectives. 


Name: 

Scan Type: 

Objectives: 

TMAP 

Use several radial scans at various 
azimuths; azimuthal scans included on 
opportunity basis; 15 cm' 1 

Optimize observations to 
divulge both radial structure 
and cooling curve 
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information. 

TDIFF 

Radial scans 15 cm 1 . Planned in pairs, 
carefully designed so unlit and lit sides 
are observed at the same geometry. 

Constrain the thermal budget 
due to vertical transport of 
heat through the rings. 

POLARIZ 

Stare at each main ring using 3 different 
orientations, with copious deep space 
interleaved. 15 cm 1 

Determine whether any 
polarization signal is emitted 
from the rings. 

EQLBN 

Radial scans 15 cm 1 

Establish a “long baseline” 
record of cooling and heating 
in the months before and 
after equinox. 

EQSBN 

Radial scans 15 cm 1 

Establish a “short baseline” 
record of heating and cooling 
in the days around equinox. 


Table 9: summary of CIRS XXM mission rings observations adn objectives 


Name: 

Scan Type: 

Objectives: 

[N,S]PXXLXX 

Radial scans at 15 cm 1 . Azimuthal scans 
on opportunity basis. 

Complete the phase/lat 
coverage started during main 
mission, extend to missed 
solar geometries (PIE’s); take 
over the previous TMAP 
objectives (non-PIE’s). 

SHADUN, 

SHADLIT 

Shadow ingress and egress scans and/or 
radial scans at 15cm ’, high latitude. 

Map the temperature 
variation around the shadow 
region for determining 
thermal inertia. 


Table 10: historical summary of CIRS rings naming conventions 


Science goal: 

Main Mission: 

XM: 

XXM: 

Composition 

COMP 




SUBM 



Radial 

SUBM 

TMAP 

[NS]PXXLXX 

Temperature 

TEMP 
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Azimuthal 

SHAD 

TDIFF 

TDIFF 

Temperature 
and Vertical 

VERT 

VERT 

SHADUN 

Transport 



SHADLIT 

Temperature 


EQSB[N,S] 


response 
around equinox 


EQLB[N,S] 


Seeking 

polarization 


POLARIZ 


Faint ring 

FMOVIE 



Observations 

FMONITOR 



Star 

occultations 

CWLEOOCC 



A one-off 
discovery mode 

ZEROPHAS 




2.3.3 Titan 

CIRS achieves different science goals at different distances from Titan (see Figure 6). 
Typically, CIRS makes the following requests (symmetric about closest approach): 

Table 11: CIRS Titan observation descriptions 

TIME FROM C/A NAME DETAILS 


o to 10 mins H IRES 

10 to 45 mins FIRLMBT 


45 to 75 mins FIRLMBAER 


Surface mapping (e.g. slew over south pole). 

1-2 radial limb scans with FPi to derive 
temperatures in the 8--100 mbar region via the N r 
N 2 collision-induced absorption between 20--100 
cm'. 

[resol: 15.5 cm 1 ; max rate: 40 prad/s; angle: -25 
mrad or - 0-250 km altitude] 

1-2 radial limb scans with FPi to measure/ 
characterize particulate and condensate 
distributions, abundances and properties. 
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[resol: 15.5 cm' 1 ; max rate: 50 prad/s; angle: -40 
mrad or - -100-700 km altitude] 

75 to 135 mins FIRLMBINT Integrate at two altitudes on the limb (125 km, 225 

km) with FPi to search for signals of CO, Fl 2 0 and 
new species, [resol: 0.5 cm' 1 ] 

2:25 hrs to 5 hrs FIRNADMAP Slow scan spanning the disk to sound tropospheric 

temperatures at 40 — 200 mbar, via the N 2 -N 2 
absorption, and to retrieve surface temperatures 
in the presence of aerosols around 520 cm *. 

[resol: 15.5 cm' 1 ; max rate: 12 prad/s, set by 1 FOV in 
5 mins] 

5 to 9 hrs (type 1) MIRLMBMAP Map 1/4 limb using the FP3 and FP4 arrays, to infer 

stratospheric temperatures via the 1304 cm' 1 band 
of CH 4 . The arrays are placed perpendicular to the 
limb at two altitudes (-350 km, 100 km), chosen to 
provide overlapping coverage of the altitude 
range 150 to 420 km. The arrays are used in blink 
(ODD-EVEN) mode. After mapping both altitudes, 
the arrays are stepped 5 degrees in latitudes for 
the next step, [resol: 15.5 cm' 1 ] 

5 to 9 hrs (type 2) MIRLMBINT As in MIRLMBMAP, except that only a single 

latitude is covered, at two over-lapping altitudes 
(-350 km, 100 km) for 2 hrs in each position. To 
search for and measure new species in the mid-IR: 
methyl, benzene etc. [resol: 0.5 cm' 1 ] 

5 to 9 hrs (type MIRLMPAIR Arrays are placed parallel to limb at a specific 

3) altitude, continuous integration for 4 hrs using 

mid-IR PAIR detector mode. Purpose is to achieve 
highest possible S/N for trace species by co-adding 
all detectors (same altitude), [resol: 0.5 cm' 1 ] 

9 to 13 hrs FIRNADCMP Integrate on the disk at emission angle 

approximately 60 degrees with the FPi detector, 
in order to measure spatial abundance distribution 
of weak species and search for new species in the 
far-IR. Use some nodding on/off disk to get deep 
space for reference/cal. [resol: 0.5 cm' 1 ] 

13 to 22 hrs MIDIRTMAP Scan the entire visible disk with the FP3/FP4 arrays 

perpendicular to the scan direction ('push- 
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broom'), to measure stratospheric temperatures 
via the CH 4 v 4 band. Used for later dynamical 
analysis, for winds, waves etc. [resol: 0.5 cm' 1 ; rate: 
max 10 prad/s; TOST type M2 includes 30-min ISS 
dwells at start and end for cloud search; M3 has 1 
ISS dwell at start or end; M4 is all CIRS] 

22 to 48 hrs COMPMAP or Map a meridian across the planet either E-W or N- 

TEMPMAP S, using the FP3/FP4 arrays in one to five positions 

longwise. To search for new species, and/or 
monitor temperatures. 

Away from TEA (Titan Long-duration (12-36 hours) stare at Titan at long 

closest approach Exploration at range. Typically one FP3/FP4 footprint. May 

Apoapse) include an ISS NAC mosaic if Titan’s angular 

diameter exceeds 6 mrad. 


25 




Mid-IR Nadir Integ 
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(composition— 

0.5 cm- 1 ) 


Mid-IR Temperature Limb Map (15.5 cm- 1 ) 


Far-IR Nadir Map 
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Far-IR Temperature Limb 
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High-res. nadir map 
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Figure 6: schematic summary of CIRS Titan observations (not all shown) 


2.3.4 Icy Satellites 

During the Cassini prime mission, there were eight flybys of the classical icy satellites 
targeted at 1000 km or less, as well as a number of 'Voyager-class' (less than 300,000 km) 
encounters. There were also flybys of several much smaller satellites, such as Janus and 
Epimetheus, at various distances. The dimensions of these objects range from less than 100 
km to as much as 1530 km (for Rhea). Consequently, it is most useful to discuss CIRS 
observations in terms of the angular diameter of the object, rather than its distance from 
the spacecraft. 

Normally the spacecraft orientation is controlled by momentum wheels (RWA), which 
provide pointing accuracy and precision of ~2 mrad and -0.04 mrad, respectively. 
Consequently, the full spatial resolution of the FP3 and FP4 pixels (0.3 mrad) cannot be 
utilized with reasonable confidence until the target exceeds 1 mrad in diameter. At that 
point, the separation of focal planes 3 and 4 (0.9 mrad), and the large size of FPi (nominally 
3.9 mrad) offer reasonable assurance of obtaining useful data in either FPi or at least one of 
FP3 and FP4. As with the planet and the rings, the order of magnitude difference in FOV 
scale between FPi and the other focal planes plays heavily in the design of the icy satellite 
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observations. Approaching from a distance, CIRS observations proceeded roughly in the 
following order. As the observation plan matured, the naming convention gradually changed 
to reflect the multi-purpose nature of the data. 

Table 12: CIRS icy satellite observation descriptions 


NAME 

FPlINT, 

FP34INT 

FPlGLOBAL 

FP34GLOBAL 

FP13STARE 


REQUIREMENTS DESCRIPTION 

1-3 cm 1 resolution Compositional integrations (typically performed at 

Q = 1-3 mrad ranges where the target angular diameter (Q) is 1-3 

mrad). Because spectral features of solid materials 
are generally broader than those of molecular 
lines, these observations are made at 1-3 cm 1 
resolution. Conducted in staring mode, using FPi 
and center detectors of FP3 and FP4, these are 
concentrated in geometries with low phase angles, 
so as to view the icy satellite surfaces at their 
warmest and maximize the signal-to-noise ratio. 
The steep rolloff of intensity with increasing 
wavenumber requires increasing integration time 
to extend the spectrum. For a 100 K surface and a 
spectral resolution of 3 cm' 1 , a SNR of 10 at 600 cm 1 
can be obtained with 1 minute of integration, but 
extending the spectrum to 800 cm 1 with the same 
SNR requires integration of 100 minutes. 
Therefore, extending the spectrum to the highest 
wavenumbers will involve coadding spectra over 
the entire tour. 


SECLN, SECLX Satellite eclipse entry (N) or exit (X) observations 

(£2>i mrad, phase angle 4> <100°, spectral resolution 
delv = 15.5 cm 1 ; a typical eclipse duration is -2 
hours). The high sensitivity of FPi permits 
observations of objects even when the FOV is 
incompletely filled, or when the phase angle is 
moderately high. At low phase angles, even FP3 
can be used to follow the initial portion of the 
cooling curve for relatively dark objects; this allows 
observations with body-centric spatial resolution 
better than 100 to be made when the apparent 
target size exceeds -3 mrad. Eclipses provide the 
thermal inertia of the upper mm or so of the 
surface, as well as estimates of surface coverage 
by relatively large fragments of consolidated 
material. 
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FPlFAZ 

FPlMAP 


FP34FAZ, 

FP34MAP, 

FP34REG 


FPGREEN 


FPlPOLE 

FPlMAP 


FPlDAYMAP, 

FPlDRKMAP 


Phase/longitude (diurnal cycle) coverage (Q>i 
mrad, Av=i5.5 cm 1 , 10-20 minutes). Focal plane 1 
will be utilized to determine the disk-averaged 
temperature of the satellites. From the resulting 
diurnal behavior, mean thermal inertias in the 
upper cm or so of the surface will be derived. 

Global thermal inertia mapping and/or hot spot 
monitoring (Q>3 mrad, Av=i5.5 cm 1 ). Maps are 
made by slewing and rastering FP3 and FP4 across 
the disk at rates not to exceed that for Nyquist 
sampling (16 prad/s in blink mode); observation 
durations will typically be 10-30 minutes. These 
maps will be successful for varying portions of a 
satellite, depending on its albedo and thermal 
inertia. The exception is Enceladus, which is so cold 
that even the subsolar regions will be barely 
detectable in an individual measurement. In this 
case, mapping will serve to monitor the satellite for 
ongoing endogenic activity; for example, active 
sources at or above the N H 3 . Pl 2 0 eutectic 
temperature will be easily observed if they fill more 
than a few percent of an FP3 or FP4 pixel. 

Search for a solid state greenhouse (Q>io mrad, cj> 
<40°, Av =15.5 cm 1 , 15 minutes). The decreasing 
absorption coefficient of water ice with decreasing 
wavenumber in the far-IR permits detection of 
radiation from increasingly far below the surface, 
reaching as deep as 1 cm at 10 cm 1 . Slow east-west 
slews using FPi allow following the daily 
penetration of the thermal wave into the regolith. 

Polar night (annual cycle) coverage (f2>io mrad, Av 
=15.5 cm 1 , 15 minutes). FPi observations of the 
dark winter polar region from high latitude permit 
determination of the seasonal cooling curve. This 
enables an estimate of thermal inertia in the upper 
tens of cm of the polar regolith. 

Flemispheric FPi mapping near closest approach 
(Q>io mrad, Av =15.5 cm 1 ). The high sensitivity of 
FPi permits rapid mapping of satellite disks near 
closest approach, where time is at a premium. 
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Nyquist-sampled maps at low spectral resolution 
(slew rate 400 prad/sec) will typically take 10-30 
minutes. These permit identification of minor 
thermal anomalies, even on the night hemisphere. 


2.3.5 Other TWTs and WGs 

The cross-discipline TWT (XD TWT) allocated time near apoapse on some revs. The 
magnetospheres TWT allocated critical time for magnetospheric observations, such as time 
near critical transition boundaries. 

2.3.6 Primes and Riders 

CIRS may operate either as the 'prime' instrument, in which case it controls the 
spacecraft attitude by a set of uplinked pointing commands, or else as a 'ride-along' partner 
when another team is 'prime'. As the mission progressed, the approach to observing also 
matured to include “collaborative” observations in which either CIRS or another instrument 
may serve as the PRIME team responsible for the developing the designs but all rider 
instruments’ pointing requirements were accommodated. Late in the mission, PIES were 
established in the science integration part of the uplink process - these observations are 
considered of highest priority in the acquisition of the data to provide some level of 
protection against deletion. The data is collected and processed in the same manner as all 
other data. 

2.4 Data Gaps 

Gaps in the calibrated CIRS data fall into two categories: data that never reached the 
ground and data that could not be calibrated. 

2.4.1 Data that did not reach the ground 

There are four primary reasons for missing data from CIRS observations: Spacecraft 
anomaly, instrument anomaly, the Deep Space Network (DSN) being unable to receive the 
data, and onboard data policing due to CIRS exceeding its allowed data volume. 

A spacecraft anomaly can present three ways. The anomaly can force the instrument to 
turn off, thus no data is taken; and an anomaly could require the spacecraft to alter its 
pointing, thus CIRS may have produced data but not at the desired target; and an anomaly 
can cause transmission of data to fail, thus preventing it from reaching the ground. 

Similarly, instrument anomalies can cause data corruption, or a complete lack of data 
output. The anomalies can be of hardware or software in nature. A hardware anomaly can 
sometimes correct itself or be corrected by the flight software, but when that does not 
occur, or when the anomaly is in the software, the solution is typically to reboot the 
instrument. 

DSN issues can be a result of hardware malfunctions with the antenna arrays, weather 
(rain and other atmospheric conditions can cause packets to be dropped), and packets are 
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typically lost at higher rates when the spacecraft is within a few degrees of solar 
conjunction. 

Data policing is rare. The standard CIRS data rate requested is 4000 bps for noadd and 
2000 bps for coadded data. CIRS flight software, however, has a data compression 
algorithm that makes the normal output between 2800 bps and 3600 bps for noadd data 
and between 1700 and 1900 for coadded data. The signal from a specific target is one of the 
primary factors in determining where in those ranges the actual data would fall. On some 
observations CIRS has concluded it would have extra data volume available after an 
observation that could be applied to a future observation. In those situations, CIRS would 
run at a higher data rate than requested and would therefore risk being data policed. On 
most occasions of data policing, it is a deep space observation run at the no-add rate but 
with a request of only 3000 bps that is curtailed. This is designed intentionally so that CIRS 
could take data for as long as possible and obtain as much data as was allowed. 

2.4.2 Data that could not be calibrated 

There are several factors that may lead to data that cannot be calibrated. Chief among 
these is interference. Interference in the raw data has been known to occur when there is 
vibration induced in the instrument scan mechanism. The most common source of vibration 
has been the spacecraft reaction wheels operating at certain angular velocities. Interference 
has also been observed during thruster firing and from malfunctioning actuators on other 
instruments (MIMI LEMMS). 

Instrument temperatures also affect the ability to calibrate the data. CIRS attempts to 
stabilize its focal plane assembly temperature; however, during heating, the temperature 
will rise and can be at a temperature that doesn’t correspond to a set-point. In those 
situations, the data is bypassed and does not get calibrated. 

Unusual laser modes can also prevent data from being calibrated. When the CIRS laser is 
turned off (generally only when the instrument is powered off or rebooted for anomalies), 
the laser can switch to a different configuration. It is impossible to know ahead of time what 
the configuration will be, but it often takes several days or weeks to stabilize. When a known 
laser mode is in use, data can be calibrated normally. However, in unknown modes or when 
the laser is drifting between modes, proper wavenumber correction of the data is 
impossible. 

2.4.3 List of Anomalies 

Below is a list of time periods and observations with significant or total data loss due to 
one of the above causes. Times listed for anomalies are the times of the last CIRS packet 
that reached the ground prior to the anomaly and the first packet received after the 
anomaly recovery had begun. 


Spacecraft Anomalies: 
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9/11/2007 5:48:35 

9/15/2007 16:30:47 

Spacecraft safing due to TWTA anomaly 

11/2/2010 21:28:49 

11/15/2010 19:30:10 

Spacecraft safing due to CDS anomaly 

Instrument Anomalies: 


10/25/2004 17:28:44 
overflow 

10/28/2004 2:39:05 

Commanding error triggered FSW buffer 

12/10/2006 4:51:21 

12/11/2006 9:18:49 

BIU Anomaly 

12/16/2006 5:43:53 

12/17/2006 9:43:20 

BIU Anomaly 

4/4/20077:30:19 

uplink 

4/5/2007 5:04:10 

Instrument turned off by accidental contingency 

5/19/2007 23:37:32 

5/20/2007 23:01:33 

BIU Anomaly 

11/2/2008 9:36:02 

11/4/2008 17:33:58 

BIU Anomaly 

2/6/2009 12:41:48 

2/10/2009 6:57:13 

BIU Anomaly 


DSN and Data Policing issues (Data Policing Issues are extremely rare on science 
observations. But since there is no simple way to distinguish them from DSN issues, they are 
lumped together here). Only observations with less than 10% of the data making it to the 
ground are listed below. 

For list of major gaps, see Appendix B. Other gaps in calibrated data may occur due to: 
Instrument interference: 

Temperature Drifts: 

Laser mode issues: 


2.5 Electrical interferences (spikes) and other data issues 
2.5.1 Background 

There are several known interference patterns in the CIRS science data. Some of them 
have been identified as relating to externally controlled hardware, others are internally 
controlled, and the source of one has never been identified. 
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For interferences external to the instrument, and regulated elsewhere on the spacecraft, 
two patterns have been identified. Both are tied to communications with the CIRS Bus 
Interface Unit (BIU), which shares a power source with the instrument and provides the only 
means of communication between the instrument and the spacecraft. 

The first interference pattern is tied to the telemetry packet collection cycle. During the 
first Instrument Check-Out test (ICOi) performed after launch, it was noticed that a pattern 
of large spikes at quasi-regular intervals was evident on the FPi interferograms. These 
occurred with intervals of 11, 11, 10, 11, 11, 10 (repeated) RTIs and were soon determined to 
originate in the Bus Interface Unit (BIU). Querying of the BIU by the spacecraft Command 
and Data System (CDS) at a frequency of 0.75 Hz to transfer cached science packet data, led 
to interferences at 1.333 second intervals, which was subsequently discretized to the nearest 
1/8 of a second (1 RTI). The CIRS telemetry rates were later redefined to be consistently 16 
RTIs apart, so from the Cooler Cover Release Checkout onwards, the interference rate was 
changed to 0.5 Hz. 

The second noise pattern that is controlled externally is tied to the spacecraft 
communications cycle; it has always been fixed at 8 Hz. Every 125 ms, the spacecraft notifies 
the instrument that is it time to check the BIU for transactions such as new commands, 
updated clock time, and to upload any housekeeping or science-engineering data. Each time 
there are transactions BIU power consumption increases linearly as data throughput 
increases. Since the science packets transactions are approximately 5 times the size of any 
other transactions, they produce greater and more sustained power transients than the 8 Hz 
transactions. As a result, the 0.5 Hz noise pattern creates stronger spikes in the science data 
than the 8 Hz pattern of the rest of the transactions. In the figures below, you can see how 
average data rates affect power consumption, and how transmissions affect current and 
voltage. 



Average Data Rate (Kbps) 


Figure 7 - BIU average power consumption. 
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SAMPLE BIU POWER 
TRANSIENT WAVEFORMS 


a) 


b) 


In-rush current from host supply 

X: 20psec/div 
Y: 200ma/div 


Vcc transients at host 
supply output 

X: 20(xsec/ div 
Y: 20mv/div 



Figure 8 - sample BIU transient waveforms. 

Complicating the issues for 8 Hz noise is an additional source of interference. This, 
second source of 8 Hz noise is one of two known noise patterns that are controlled by 
timing established in the CIRS flight software (FSW). The first FSW regulated noise follows 
the schedule of the numerical filtering and compression cycles of the instrument. 
Traditionally, this timing was linked by the FSW to the same 8 Hz cycle used for 
communications. Through on board experimentation, it was discovered that this internal 8 
Hz noise could be separated from the externally controlled 8 Hz noise. The FSW was 
updated in July 2010 to implement that separation so that the internally controlled noise 
could be averaged out over multiple interferograms. 

The other source of noise with a FSW controllable rate is noise caused by performing 
analog MUX data readouts. The noise was discovered after a FSW change caused internal 
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timing to be more tightly controlled. As a result, the noise no longer averaged out over a 
large number of scans, it appeared in the same positions on each scan. The resulting fixed 
pattern produced a spike in the spectra. This noise is the weakest of the known noise 
patterns and is nearly impossible to see in uncalibrated instrument data. 


Figure 9 below is a plot of averaged interferogram data from focal plane 1 in early 2005. 
The 0.5 Hz spikes stand out on the sides, with the fifteen 8 Hz spikes between them (The 
sixteenth 8 Hz spike coincides with the 0.5 Hz spike so it is masked.) The MUX related noise, 
which would have been observed in individual interferograms at 0.77-1.0 Hz is not 
discernible. 


Averaged FP l_Interferogram_Data 

2005-032T00:28:06 - 2005-032T05:58:06; Sample Count: 2764; scanjength: 400; Shutter: OPEN; 

-15.0 ; : 

-17.5 

-20.0 

-22.5 

-25.0 

-27.5 

-30.0 

-32,5 

-35.0 

-37.5 



-47.5 S y 

-50.0 

-52.5 

-55.0 

I 

-57.5 

-60.0 

-62.5 

-65.0 

-67.5 

-70.0 

-72.5 

-75.0 1 •- - - 1 — -j — 1 - - -4 i | | ■ 1 - 1 - - \ 

1,600 1,610 1,620 1,630 1,640 1,650 1,660 1,670 1,680 1,690 1,700 1,710 1,720 1,730 1,740 1,750 1,760 1,770 1,780 1,790 1,80C 

Sample Number 

| — Interferogram Data] 


Figure 9 - plot of averaged FPi interferogram data showing electrical interferences (spikes) 


Each of the above noises has been observed in focal planes i and 3. The internal filtering 
related noise has also been observed in FP4 during a noise experiment, but no noise has 
been seen in FP4 under normal operating conditions. 
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A further ‘spike’ pattern only affects FP3. It occurs at a rate of approximately 8.33 Hz. 
The CIRS team has performed experiments in attempts to affect the timing of the noise; 
however, none of the experiments have had any affect. All that can be learned was that the 
timing of the noise is relative to the start of the scan. 

In addition to spikes occurring in the spectra at positions corresponding to these 
frequencies, additional spikes occur at the harmonics of each frequency. For example, in FPi, 
the 8 Hz noise appears at 191.38 cm 1 , but it also appears at 382.62 cm 1 and 573.93 cm' 1 , 
corresponding to the 16 Hz and 24 Hz harmonics. While easily identifiable, the spikes in the 
spectra are not the same size for each harmonic. The conversion between the frequency of 
the noise in the interferogram and the location in the spectra of the corresponding spike is 
approximately 23.91375 cm 1 per 1 Hz (though it is affected by the laser mode which was 
active at the time the data was recorded). 

The final interference pattern does not consist of a spike comb. Two sine waves (at least) 
are present: one in FPi and one in FP3 interferogram data. This interference presents itself 
as large spikes in the spectra. Its frequency changes over time, so the spike moves within 
the spectra. For FPi it has been seen in the range -23 - 430 cm 1 range (at least), and for FP3 
between 880 and 1020 cm 1 (at least). The FP3 sine wave spike is a harmonic of the one in 
FPi. The exact cause of the sine waves is unknown; however, it is known that the 
frequencies decrease as the CIRS electronics board temperatures increase. 

2.5.2 Mitigation 

Over the years, several steps were taken to eliminate or average out the noise spikes. In 
2005, CIRS began using scan lengths that would create 16 positions for the 0.5 Hz noise 
spikes (39, 57, 97, 225, 8( 401 RTIs). Algorithms were also introduced into the ground data 
processing to (i) locate the spikes, and then (ii) to subtract out a typical shape (‘comb’) of 
the noise from the interferogram, thereby greatly reducing the appearance of the noise. 

In FSW Version 6.0.0 (uploaded July 2010), the ability to manipulate the timing of events 
that affect the position of noise spikes relative to the start of an interferogram was 
introduced. To that end, the noise spike that occurs in FSW Version 5 at a 1 Hz rate, and the 
stronger of the two observed {8, 16, 24} Hz interference sources, have their rates 
randomized. The positions of the noise spikes that occur at 0.5 Hz, as well as the positions 
of a second, weaker, noise pattern occurring at (8, 16, 24} Hz, are manipulated by adding 5 
ms to the scan length. For example, in versions of the FSW previous to FSW 6.0.0, a 401 RTI 
scan was approximately 50.125 seconds in duration; in FSW 6.0.0 the duration is 50.130 
seconds. 

The result of this additional time being added to each scan is that the 0.5 Hz and {8, 16, 
24} Hz noise patterns are shifted within the interferogram relative to the start of fly back, 
and thus to the start of scanning. The 5 ms shift quantity creates 200 potential time delays 
from the start of the scan to the first 0.5 Hz noise spike in the interferogram. Assuming 
perfect timing, when 200 consecutive scans are averaged the combined noise pattern would 
appear to be a 100 Hz pattern. With the {8, 16, 24} Hz noise, there would be 25 potential 
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positions inside the interferogram for the noise, which would produce a 200 Hz noise 
pattern when averaged. Both the 100 Hz and 200 Hz noise patterns lie outside the 
bandwidth of CIRS Focal Planes 1 and 3 and therefore do not show up in the interferogram. 
The full amplitude 0.5 Hz spikes and a much diminished {8, 16, 24} Hz spike sequence are 
present in each individual interferogram. 

In summary, the new timing configuration was as follows: 

• Science packet collection, 'A Hz noise, shifted 5ms per scan 

o Individual scans: 200 noise positions; but Vi Hz noises still present in each 
individual scan. 

o Averaged scans: 200 positions average to a 100 Hz noise; outside the band pass 

• Analog MUX channel data collection, 1 Hz noise 

o Individual scans: randomized noise positions; noise outside band pass 
o Averaged scans: randomized noise positions; noise outside band pass 

• Spacecraft communications cycle, 8 Hz noise 1, shifted 5 ms per scan 

o Individual scans: 25 noise positions; but 8 Hz noise still present in each individual 
scan. 

o Averaged scans: 25 positions average to a 200 Hz noise; outside the band pass. 

• Internal/numerical filter related, 8 Hz noise 2. 

o Individual scans: randomized noise positions; noise outside band pass, 
o Averaged scans: randomized noise positions; noise outside band pass. 

• lnternal/FP3 only, -8.33 Hz noise 

o No change - still always 8.33 Hz 

The effects of the changes were positive and are shown below. Averaging scans with 
noise in the above configurations reduces the spikes to levels smaller than the already 
despiked data from earlier versions of flight software. 


36 



5 

Total IFGMs in FSW 5.0.5 block: 7637; Total IFGMs in FSW 6.0.0 block: 10484 


FSW 6.0.0 - FP1 Averaged FIFM (9607 scans; August 1 - September 22, 2010) 


FSW 6.0.0 -FP1 Averaged FIFM (6098 scans; NPTS 2775-2785; August 1 - September 22,2010) 


FSW 6.0.0 - FP1 Averaged IFGM (6098 scans; NPTS 2775-2785; August 1 - September 22, 2010) 


FSW 5.0.5 - FP1 Averaged FIFM (6331 scans; May 1 - June 1, 2010) 


— • FSW 5.0.5 -FPl Averaged FIFM (6098 scans; NPTS 2775-2785; May 1 - June 1,2010) 


FSW 5.0.5 - FPl Averaged IFGM (6239 scans; NPTS 2775-2785; May 1 - June 1, 2010) 
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Figure 10 - comparison of averaged FP4 interferogram data from different mission periods 
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In the FPi and FP3 average spectra the results were even more pronounced. Note that the 
large 8.0 Hz spikes in the blue line (pre-FSW 6.0) are not present in the red line (post FSW 
6.0): 


Average Power Spectrum - FPI 

209)Q-ffi28T0Q:Q4:OQ - 2O!t0-(D28TO8:lS:O8; Sample Count: 2778; scan_len_calc: 401; Shutter: OPEN; 



ower Spectrum | 


Figure 11 - FPi uncalibrated s pectrum before (blue) and after (red) despiking 
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Similarly in FP3, the 8 Hz spikes are essentially gone, while the 8.33 Hz spikes (598, 797, 
996) are visible, as expected, in data acquired with both versions of the flight software. 

Average Power Spectrum - FP3_5 



ower Spectrum | 


Figure 12 - raw FP3 uncalibrated spectrum before (blue) and after (red) despiking 

The solution is not, however, perfect. By randomizing the positions of some spikes and 
shifting others within the interferograms, the standard deviations of the samples at each 
position are increased. In Figure 13 below, you can see that the 8 Hz spikes in green (pre- 
FSW 6.0) are much larger than the 8 Hz spikes in black (from a test of the FSW 6 
configuration). A similar statement can be made about the pre-FSW 6.0 standard deviation 
in red and the standard deviation of the FSW 6.0 test, in purple. However, the baseline of 
the new configuration’s standard deviation is higher in the regions without a spike. 
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Figure 13 - Saturn spectrum before and after despiking 

Note: since both the amplitude and phase of the noise spikes are variable, it is usually 
necessary to confirm the presence of an apparent artifact in the spectrum by looking for a 
spike in the standard deviation at the same wavenumber position. The effect is especially 
clear on averages of deep space spectra. 
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3-0 Calibration of the instrument 


3.1 Overview of data flow 

The block diagram in Figure 14 schematically illustrates the overall flow of data from the 
spacecraft, through various transfers on the ground, to reach the user. In this document we 
will focus only on the part of the data pipeline that occurs at GSFC (top center dashed box). 



Figure 14: block diagram of data flow from the CIRS instrument to the eventual user. 


A more detailed picture of the data flow at NASA GSFC is shown in Figure 15. 

Raw compressed instrument data packets arrive from the spacecraft, via the DSN, at 
NASA J PL. An automated query program brings all new or changed packets via secure copy 
to the GSFC SOPC computer, shown in the upper left behind the firewall. At this point they 
are requeried to an archive computer, and across the firewall to the main data server 
computer (cirsdata). Decompression of the data to form raw table blocks is performed here 
automatically every night. Following the creation of the raw interferograms (IFM), pointing 
and geometry data is automatically generated based on spacecraft ephemeris data (see 
Section...). At this point the level 0 ‘raw’ database is complete. 
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Figure 15: flowchart of CIRS data at GSFC. For definition of Vanilla see Section 4.1.2 

The second part of the data flow takes the level o (integer count) interferogram data 
and filters out bad data, applies some ‘cleaning’ algorithms to the good data, creates a 
record of any potential issues for calibration (e.g. scan mechanism out of lock), and re-writes 
the interferograms as floating point data. See Section 3.3 for further details. 

The third and final process is the calibration, which combines interferograms of 
reference temperature sources (space for FPi, space and shutter for FP3/4) with the on- 
target data, and then calibrates both the radiance (from counts to physical units) and the 
wavenumber axis (from sample number to inverse centimeters). See Section 3.4 for details. 

3.2 Instrument readout to raw tables 

CIRS data is normally read out at a rate of 4 kbits per second (kbps) 4 and 8 kbit packets 
are transferred via the Bus Interface Unit (BIU) to the spacecraft Command Data Subsystem 
(CDS) every 2 seconds. CIRS has many ways of reducing data rate, such as co-adding IFMs 
(reduces rate by factor 2), dropping FP3 or FP4 readout (reduces by 45%), dropping FP3 and 


4 6 kbps during the Jupiter encounter. 
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FP4 readout (reduces by 91%), or going to housekeeping only mode. Data packets are 
compressed in the instrument electronics, packetized on the spacecraft, transmitted to the 
Earth via the Deep Space Network (DSN) and then decompressed by software on the 
ground. At this stage, the data goes into a processing pipeline, which organizes the science 
data into binary tables and interpolates housekeeping records. Pointing and geometry 
tables are also produced when the NAIF data becomes available. 


3.2.1 Raw data tables 

The raw data tables are fixed or variable-length binary tables in the Vanilla format, 
created by simply uncompressing the instrument packets from the spacecraft and writing 
the data as binary files. There are four types: 

• OBS : observation data. A fixed-length record binary table containing basic 
information about the data modes during the observation, including which 
detectors are active, whether COADD mode is activated, the scan length (spectral 
resolution) and so on. 

• IFGM : variable-length binary records containing the raw interferogram data, in 
units of detector counts. 

• FRY : fringe voltage data. 

• IHSK : interpolated housekeeping data, including the temperatures recorded by 
sensors on various parts of the instrument, which is needed for eventual 
calibration, housekeeping data packets are created every 64 seconds, however to 
align this information with the individual interferograms (4.5 to 52 seconds in 
length), the housekeeping data is interpolated to match the start times of each 
interferogram. 

Further information about the Vanilla tables, how to read them and their contents, is 
given in “Section 4.0 How to search and extract data”. 

3.2.2 Navigation data tables 

The NAV tables take as input various kernels provided by the Cassini Navigation Team 
and produce various quantities related to the location and velocity of Cassini and other 
bodies in the target system and the orientation of CIRS. This is done almost entirely by using 
the NAIF toolkit, a collection of software which is written by JPL’s Navigation and Ancillary 
Information Facility, which reads the various formats that the kernels come in and performs 
various operations on that data. (A list of the kernels used for every NAV table and the 
version of the NAIF toolkit that was used to produce a given table is entered in the text .LBL 
file associated with each binary .DAT file which contains the actual table.) Important types 
of kernels are the SP kernels, binary kernels which contain the ephemeris for Cassini, the 
Earth, the Sun, and the bodies in the target (either Jupiter or Saturn) system, the C kernels, 
binary kernels which contain information about the orientation of Cassini at any given time, 
and the PC kernels, text kernels which contain the masses, rotation rates and axes, shapes 
(as modeled by triaxial ellipsoids), etc., for the bodies in the target system. Other kernels 
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contain the details of the spacecraft and instrument coordinate systems, and the 
relationship between spacecraft event time and the “spacecraft clock”. 

It is important to keep in mind that the software that produces the NAV tables has no 
knowledge of what a particular observation is for or what the intended target is. For that 
reason, entries in the NAV tables may seem to be extraneous. For example, an entry in the 
RIN (rings) table is always made, even if the rings aren’t a target and even if CIRS is pointed 
away from the rings. 

Another important point to keep in mind is the use of NAIF reference ellipsoids by the 
NAIF toolkit. For bodies with solid surfaces, those surfaces are modeled by a triaxial 
ellipsoid, with axes contained in the planetary constants (PC) kernels. For gas giants (Saturn 
and Jupiter), the NAIF reference ellipsoid is suppose to model the i bar surface. Flowever, 
using the reference ellipsoid as the one bar surface is a very poor choice, since, because of 
the subtleties of the atmospheres (winds, deviations from axisymmetry, etc.) the NAIF 
reference ellipsoid can deviate from the actual one bar surface by over 100 km. The NAV 
tables always use the NAIF reference ellipsoids as references, but keep in mind they may be 
inadequate for the problem you are trying to solve. 

The CEO (geometry) tables 

The GEO tables are unique among the NAV tables in that they do not depend at all on the 
orientation of Cassini, but only upon its state vector (location + velocity) and the state 
vectors of the Sun and all of the bodies specified in table 4.4.1. The list of quantities included 
in the GEO tables can be found in the file GEO.FMT, and are more fully described in 
DATASIS.CAT. The field keys for the GEO tables are scet and bodyjd. bodyj'd is similar to 
targetjd for the POI, RIN, and TAR tables, but is called bodyjd because there is no “target” 
to speak of. All of the bodies in table 4.4.1 have an entry in the GEO tables for every scet. 
The entries in the GEO tables are independent of a particular CIRS detector. 

Note that the quantities in the GEO tables are computed not at the scet, but at the 
time_zpd, the time of Zero Path Difference where the majority of the signal is collected, 
typically a few seconds after the scet. Also note that while the .LBL files for the GEO tables 
list C kernels in the collection of kernels used to create an individual table, the C kernels are 
not actually used. 

The POI (pointing) tables 

The POI tables contain quantities that depend on the orientation of CIRS and on the 
individual detectors. The three field keys are scet, det, and targetjd. (Flowever, as with the 
GEO tables, many of the quantities in the POI tables are actually computed at time_zpd, and 
a few at time_end, the end time of the scan.) The entries in the POI tables are found in the 
file POI.FMT. 
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Entries are made in the POI tables for all active detectors and for every body that is in the 
field of view of an active detector. One exception is the primary (Saturn), for which an entry 
is always made at every scet and det, whether or not it is in the field of view of a detector. 
“Field of view” also requires a little explanation. Every body gets a “penumbra” attached to 
it before it is determined whether it is in the field of view or not. For airless bodies, a 
penumbra of 0.001 mrad is added to the three axes of the triaxial ellipsoid defining the 
surface of the body. For bodies with atmospheres (Saturn and Titan), the penumbra is 
larger, and depends on the detector. For FPi, the penumbra is 600 km + 0.004 mrad (i.e., 
600 + [0.004 * distance to target] is added to each axis). For FP3 and FP4, 0.003 rnrad is 
used in addition to the 600 km. Then we check to see whether this new object (target + 
penumbra) is in the field of view. If any Q point (grid of nine points defining the outline of 
the detectors: see Figure 27 ) intersects the object, then it’s in the field of view. If the entire 
object is within the field of view but not touching any Q points, then it’s within the field of 
view. Otherwise, it’s out of the field of view. Note that, even though an entry is always 
generated in the POI table for Saturn, it is only marked as in the field of view if it actually is. 

The POI tables have many entries that are computed at each Q point for the active 
detectors, but some are computed only for the boresight of the detector (Q point 5). Two 
quantities, filling_factor and emission_angle_fov_average merit special attention because 
they are computed by doing a Monte Carlo estimate. The field of view of the detector in 
question is filled with rays randomly chosen in the field of view, and then we determine 
whether or not the ray intersects the target in question. Unlike the field of view 
computation for the Q points, we do not add the penumbra and use the NAIF reference 
ellipsoid as the target’s surface. Then the filling_fa ctor is computed by taking the ratio of 
rays that intersect to the total number of rays, and emission_angle_fov_average simply 
averages the emission angle for the rays that intersect. There’s an additional complication 
for FPi. FPi is more sensitive toward the center than toward the edges. This can be 
modeled as a Gaussian response, and we use this Gaussian to bias the selection of the 
random rays toward the more sensitive part of the detector in the center. So while 
filling_factor for FP3 and FP4 is simply an approximation to the area of the detector that is 
filled with the target, for FPi weight is given to the more sensitive part of the detector. For 
purposes of speed, a modest number of random rays is used. We stop whenever 1000 
points have intersected. If 5000 points have been tried before 1000 points have intersected, 
both filling_factor and emission_angle_fov_average are set to -200, an otherwise impossible 
value. So the formal statistical error for these quantities is on the order of a few percent. 

The smear is defined in the following manner. The position of the boresight on the target 
at time_zpd is computed and saved. Then, at time_end, the angle between that point and 
the position of the boresight at time_zpd is computed, smear is the ratio of this angle with 
the angular width of the field of view of the detector, or in other words, is the fraction of the 
field of view over which the original intersection point of the boresight has moved between 
time_zpd and time_end. If the boresight at either time_zpd or time end is off the target, 
then this is set to -200. 
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Latitudes in the POI tables are given in two systems, planetocentric and planetographic. 
For spherical bodies they will be identical. 

The POI tables are also unique in that they may contain entries for targets that are fixed 
stars, which are used occasionally for testing. The list of potential fixed star targets is in 
Table 13 . Since many of the entries have no meaning when referred to a star (for example, 
latitude and longitude), these are set to the usual “CIRS bad value”, -200. Stars were given 
artificial target jd’s which have nothing to do with actual NAI F id’s in this case. 


Table 13: a list of the possible stellar targets that could be used for testing purposes that could have entries in the POI 
tables. The numerical id's are artificial, and have nothing to do with NAIF id's. Not all of these were used during the 
mission. 


bodyjd/targetjd 

Object 

10101011 

VYCMa 

10101012 

eta Car 

10101013 

CW Leo 

10101014 

AOri 

10101015 

A Sco 

10101016 

R Dor 

10101017 

0 Ceti 

10101018 

NMLTauri 

10101019 

W Hya 

10101020 

Gamma Cru 


The R IN (rings) tables 

The RIN tables contain quantities describing the rings of Saturn as the target. Like the 
POI tables, the RIN tables depend on the orientation of CIRS aboard Cassini. Unlike the POI 
table, no effort is made to determine whether or not the rings are actually in the field of 
view of a detector. An entry in RIN is made for every active detector at every scet. The field 
variables are scet and det. As required by the Planetary Data System, the coordinate system 
used to describe the rings is an inertial coordinate system centered on Saturn, with the z axis 
aligned with Saturn’s rotational axis. The rings are assumed to be in the z=o plane, out to 
140270 km. Longitudes are measured eastwards from the ascending node of the 
“intersection of the ring plane with J2000”, or, in other words, the intersection of the ring 
plane with a plane parallel to the ecliptic plane of the J2000 coordinate system. 
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The TAR (target) tables 

The TAR tables answer a simple question, whether a given object is in the field of view of 
a detector at a given scet. The field variables are scet and det. The TAR tables consist of 
single byte entries containing either a one (true) or a zero (false) for each body in 

Table 14 or Table 15, a single byte containing either one or zero which says whether 
any of the stars in T able 1 3 are in the field of view, and a byte that says whether or not the 
rings (Saturn or Jupiter) are in the field of view. In addition, the TAR tables contain 
fov_targets, an integer in which each bit is a separate entry for each body identical with the 
byte fields. The bit values are fully described in DATASIS.CAT. 

In the TAR tables, a satellite that is being eclipsed by its primary (e.g. Titan behind 
Saturn) is not marked in the field of view. However, a satellite being eclipsed by another 
satellite will still be marked in the field of view; i.e. if Enceladus is completely behind Titan, it 
will still be marked in the field of view if it otherwise would be if Titan weren’t there. 

Table 14: list of objects that were included in the NAV tables during the Jupiter flyby. The body/target id is identical with 
the NAIF id of the body. 


body_id/target_id 

Object 

501 

lo 

502 

Europa 

503 

Ganymede 

504 

Callisto 

599 

Jupiter 


Table 15: list of objects that are included in the NAV tables from the Saturnian system. The body/target id is identical with 
the NAIF id of the body. 


body_id/target_id 

Object 

601 

Mimas 

602 

Enceladus 

603 

Tethys 

604 

Dione 

605 

Rhea 

606 

Titan 
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body_id/target_id 

Object 

607 

Hyperion 

608 

lapetus 

609 

Phoebe 

610 

Janus 

611 

Epimetheus 

612 

Helene 

613 

Telesto 

614 

Calypso 

615 

Atlas 

616 

Prometheus 

617 

Pandora 

618 

Pan 

699 

Saturn 
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3.3 Cleaning and filtering of interferograms 
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Figure 16: block diagram showing checks performed on raw interferogram data 
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3-3-1 Overview 

The first stage in the conversion of an FPi, FP3, or FP4 raw interferogram (IFGM) to a 
filtered interferogram (FIFM) suitable for calibration is to send the IFGM through the CIRS 
noise detector algorithm. This algorithm tests each IFGM for disturbances from external 
noise effects, such as transient spikes, ZPD shifts, laser mode drifts, and scan mirror velocity 
variations. The logic of the algorithm is illustrated in Figure 16 (above). 

Once an FPi I FGM passes successfully through the noise detector with a value of NOISE = 
1 (NORMAL IFGM), it is sent to the CIRS FPi spike suppression algorithm. If the IFGM is 
successfully de-spiked by the spike suppression algorithm it is written to Vanilla tables as an 
FIFM. FP3 and FP4 IFGMs are presently converted directly to FlFMs by simply converting 
their amplitudes from integer to floating point format, and they are then written to Vanilla 
tables. A summary of the issues involved in suppression of FPi spikes follows; a fuller 
description of the present spike suppression algorithm is given in Appendix C. 

3.3.2 Suppression of FPi Spikes: Introduction 

Interferograms from the CIRS FPi and FP3 detectors are contaminated by two series of 
prominent 0.5 Hz and 8.0 Hz noise spikes that are caused by electrical interference due to 
real-time interrupt (RTI) signals from the Cassini spacecraft bus interface unit (BIU). The 
spikes are spaced in time during every interferogram by the RTI period of 1/8 second, and 
their intensities are enhanced every 2 seconds by a data transfer event. The phase of the 
spike pattern in each interferogram, i.e., the position of the spikes relative to the zero path 
difference (ZPD), depends upon the starting time of the scan, which can be controlled to 
some extent through selection of the interferogram length. The 0.5 Hz spike pattern in the 
interferogram results in a series of spikes in the spectrum spaced by 12 cm' 1 , corresponding 
to the harmonics of 0.5 Hz. The 8.0 Hz spike at 191.38 cm ' also appears in FPi spectra as 
strong harmonics at 16.0 Hz (382.76 cm 1 ) and 24.0 Hz (574.14 cm '). Both the 0.5 Hz and 8.0 
Hz spikes coincide with certain spectral features and can make atmospheric retrievals 
difficult. In addition, FPi far-infrared interferograms are contaminated by a prominent single 
frequency (“sine wave”) feature of unknown origin that has been observed between 18 cm ' 
and 400 cm 1 . FP3 interferograms and spectra also exhibit a harmonic series of 8.3 Hz spikes 
(see Section 2.5 Electrical interferences (spikes) and other data issues” for further details). 

The CIRS FPi multi-parameter spike suppression algorithm utilizes high fidelity empirical 
0.5 Hz and 8.0 Hz spike waveforms (“combs”) created from real deep-space data using large 
averages of FPi interferograms. The algorithm matches the positions and intensities of the 
spike combs to the observed spikes in each individual interferogram and performs checks to 
optimize and characterize the quality of the fitting process. Before de-spiking, each 
interferogram is over-sampled by a factor of 8. 

The CIRS noise suppression algorithm has proven extremely successful in reducing 0.5 Hz 
and 8.0 Hz electrical noise spikes, as well as ripples < 50 cm 1 and > 300 cm 1 , in a variety of 
Focal Plane 1 spectra calculated from interferograms recorded at 96, 97, 224, 225, 400 and 
401 RTI scan lengths. 
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Since October 2004 (Titan observations) and March 2005 (Saturn and all other 
observations), CIRS has been employing an additional technique for suppressing 0.5 Hz 
spikes. The 0.5 Hz spike pattern is effectively suppressed in averages of > 16 scans by 
recording interferograms with 16N+1 RTI, e. g. 97, 225, and 401 RTI. The 0.5 Hz suppression is 
affected by the 0.5 Hz spike pattern being averaged into the 8.0 Hz pattern in the course of 
the averaging. Unfortunately this technique is not applicable to the 96, 224, and 400 RTI 
interferograms that were acquired before the above dates. In addition, the 8.0 Hz spikes 
remain in the spectra. 

It was necessary to modify the FPi 0.5 Hz and {8.0, 16.0, 24.0} Hz spike suppression 
algorithm after CIRS flight software (FSW) Version 6.0.0 was uploaded and implemented on 
July 28, 2010 (see Section “2.5.2 Mitigation”). FSW 6.0.0 significantly reduces 8, 16, and 24 
Hz electrical interference noise spikes in both individual and average FPi interferograms. 0.5 
Hz noise spikes are still present at their pre-FSW 6.0.0 shapes and amplitudes in individual 
interferograms, but their positions relative to the start of each scan are changed. 

In light of the changes to the spike patterns resulting from FSW 6.0.0, the FPi spike 
suppression algorithm was modified to bypass {8.0, 16.0, 24.0} Hz spike suppression for each 
NOADD interferogram recorded post-FSW 6.0.0. As with earlier FSW versions, the algorithm 
subtracts out the interferogram’s D. C. baseline and deep-space shape function prior to 
suppressing the 0.5 Hz spikes. Due to the difficulty of discerning 0.5 Hz spike patterns in 
individual and average post-FSW 6.0.0 FPi COADD interferograms, both 0.5 Hz and {8.0, 
16.0, 24.0} Hz spike suppression are bypassed for COADDs. The spike suppression algorithm 
only subtracts out each COADD interferogram’s D. C. baseline and deep-space shape 
function. 

3.4 Calibration of the data 

This section describes the wavelength and radiometric calibration of the data. 

3.4.1 Wavelength calibration 

Wavelength calibration is achieved by use of a reference diode laser (783 nm) 
interferometer that uses the same scan mechanism as the IR interferometer. The scan 
mirror moves at a constant speed of 0.0208 cm/s. While the mirror is moving, data are 
logged at every zero-crossing of the laser intensity, i.e. every 391.5 nm of light path 
difference, or every 195.75 nm of mirror travel. 

Number of samples returned 

FP3 and FP4 in the mid-infrared each sample from five detectors out of a possible ten, 
depending on the selected observation mode (see Section “2.1.3 Instrument operation 
modes”). Each detector is sampled in turn, so that the number of points per mid-IR 
spectrum is 5X less than in the far-IR (FPi) spectrum. Finally, all the channels are further 
decimated by resampling factors of 18, 4.5 and 5 respectively for FPi, FP3 and FP4, so as not 
to exceed the necessary Nyquist sampling rate for the spectrum. 

The final number of points in the spectrum is therefore: 
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light path change x 1 (l) 

y logging interval resample factor 

(2x0.0208xt) 1 , v 

n„ = (2) 

(783/2X10 ) (^d Xri resamp) 

where rid is the number of detector elements switched between, and n re samp is the 
resampling factor given above. For example, aios scan with FP4 would have ^=425. In fact, 
the above formulae only provide an approximate guide to the number of samples that are 
returned from a scan of given length. This is due primarily to mirror fly-back time (at start of 
scan) which shortens the effective scan duration, and moreover requires an amount of time 
that is variable depending on the length of the preceding scan - i.e. where the mirror was 
left prior to the current scan. In addition, loss of scan mechanism phase lock due to 
spacecraft vibrations (mainly from the reaction wheels) and inherent velocity variations in 
the scan mechanism can lead to variations of ~i% in the amount of samples returned. 


Sample point spacing 

The spacing of output points on the power spectrum is the reciprocal of the maximum 
path difference of the beam, X: 

Sv = — ; X = 2x0.0208x(t- t„) (3) 


where t 0 is the time after scan start when the mirror reaches the zero path difference (ZPD) 
point, equality in the light path in both interferometer arms. This is triggered by the white 
light fringe of a reference source. Alternatively, substituting for t in Eq. (2) we arrive at: 


Sv 


1 


2XCppX7l L 


(4) 


where n L is the actual (returned) number of samples in the long side of the interferogram 
(after ZPD), and we define the constant Cfp for each focal plane as follows: 

Cpp = Tld X?2 resam p X — X 10 (5) 

For example, Cfp 4 =9. 7875x10 4 cm' 1 , so if there were 350 actual samples after ZPD in a 10 s 
scan on FP4, the point spacing in the final power spectrum would be 1.46 cm' 1 . Note that this 
Nyquist point spacing is about half the spectral resolution, depending on definition. So in 
this example the spectral resolution is about 2.92 cm' 1 . 


In practice, Fourier transforms are often made using an array size that is some power of 
2: i.e. 512, 1024 etc. Therefore, the zero-padding that occurs to the interferogram in the array 
results in an oversampling on the resultant power spectrum. In this case, the 350-point array 
would be padded to 512 and the resulting point spacing would be 0.998 cm' 1 , less than half 
the resolution. This topic is discussed in greater detail later (Section “3.5.1 ”). 

Laser modes 
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A further complexity is added to the wavelength calibration due to the fact that the 
diode laser may ‘mode-hop’ to different frequencies. This occurs infrequently, and modes 
generally persist for months or years. However even within the modes there is some smaller 
drift over time that may be due to thermal effects. Modes are generally determined by 
monitoring the fringe voltages returned in the housekeeping data packets. Therefore, care 
must be taken to match laser modes/fringe voltages of target scans, and blackbody 
reference scans, during calibration. See Figure 17 for an overview of the modes. 


Reference Interferometer Signal - Laser A: 2004 - 2011 @ 3 Secs 



gggggggggggggggggg 

000000000000000000 


Figure 17 : CIRS laser modes 

Further discussion of laser modes is beyond the scope of the User Guide. 

3.4.2 Radiometric calibration 

Radiometric calibration in the far-IR (FPi) is simplified due to the fact that there is only 
one temperature inside the instrument. Therefore, use of a deep space (2.7 K) reference 
target is sufficient to calibrate radiance. The fact that the mid-IR has two temperatures to 
contend with means that a second reference target is required. Therefore, in the mid-IR 
optical path, a shutter may be lowered for calibration purposes, which gives a 170 K 
reference body. Due to phase differences that may arise in the beamsplitter, differencing of 
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target-space and target-shutter is conducted in the interferogram domain, prior to Fourier 
transforming. 

A somewhat idealized description of the radiometric calibration is included hereafter. In 
the real pipeline, additional complexities must be accounted for, including apodization to 
reduce ringing (Gibbs Phenomenon) in the final result. 

3.4. 2.1 Idealized Interferometer Model 

Consider the following simplified radiometric model for CIRS, from ( Brasunas, 2002 ). 



| moving 
O I mirror 


the instrument 


Figure 18: simple model of interferometer 

Instrument self-emission includes emission from the detector, and from the beam 
splitter (B/S) if it is absorbing. As the optical path difference 5 (6 is the two-way difference, 
expressed as a phase) changes, the received signal (interferogram) at the detector changes. 

Let the complex reflection amplitude r of the B/S be r exp(iSr), and the transmission 
amplitude t be t exp(ibt) (assume either TM or TE polarization). Let a^ be the propagation 
amplitude from port 1 back to port 1 , thus: 

a n =( r 1 exp(/2<5r)+ r exp(i2St)exp(iS) )exp(/<5 0 ) (6) 
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where S 0 is round-trip travel common to each of the two arms. Thus the fraction of the 
input power returned to port i, P 11; is: 


P n = | ai] p=i? 2 + T 2 +2RTcos{d + 2(dt -dr)) (7) 


where R = r 2 and T = t 2 (r, R, t and T are all real). Similarly, the propagation amplitude a 12 


from port i to port 2 is 

a n = a 21 = rtexp(i(dr + <5f)Xl + ex p(/'<5))exp(z<5 0 ) (8) 

and the fraction of the power transferred from port 1 to port 2 is 

Pn = Pi, = 2RT (cos(<5)+l) (9) 

Lastly &22 an d p 22 are 

a 22 = ( r 2 exp(/2<5z')exp(/<5)+ t 2 exp(i2dt) )exp(z<5 0 ) (10) 

P 22 = |a22f = R 2 + T 2 + 2RT cos(d - 2 (dt - dr)) (11) 


Eqs. ( 6 ) through (11) are not restricted to a loss-less dielectric, but they do assume a 
specific polarization state. For a loss-less dielectric (St - Sr) is jt/2, thus the two modulated 
output signals due to the target are complementary, as can be seen from and P 12 . 
Define A to be 2(St - Sr), which is 7 t for a dielectric. 

Now let the B/S be absorbing. Consider radiant flux emitted into the FTS from port 2: the 
total fraction transmitted or returned is P 21 + P 22 - Therefore the absorption of the B/S on 
side 2, A 2 , satisfies the equation 


a 2 - 0 “ P 21 - P 22)- 


( 12 ) 


Similarly the side 1 absorption of the B/S is 

Ai=( 1 - p i 2 - p n)- 03) 

By Kirchoff’ s law, the emissivity equals the absorption coefficient. Now consider the flux 
balance at port 2: flux is lost (the flux from port 2 that is not returned) and flux is gained 
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both from port i and from the B/S emission. Let the unity emissivity fluxes from target, 
detector and B/S be Ft ar , F^et and Fg/s, a H assumed to be in thermal equilibrium and 
therefore functions of their respective temperatures (through the Planck function). The flux 
net gain or loss dF at port 2 is: 

dF = F tar P n + F bis (l - P 2l -P 22 )~ F det (\-P 22 ). (14) 

Using Eqs. ( 9 ) and ( 11 ) in Eq. ( 14 ), and considering only the modulated part dF moc j of the 
net flux (dF is a narrow spectral slice of the modulated part of the interferogram), we have 

dF moi = 2RT(F lar cos 6 - F B/S (cos 5 + cos(A - A)) + F ia cos(<5 - A)) , (15) 

dF mo d = 2 RT{(F tar - F b/s )cosS + (F det - F b/s )cos{S - a)) (16) 


Thus the self-emission anomalous phase (when A is not 7 t) becomes associated both with 
the B/S emission and with the detector emission. This model correctly shows zero net flux 
when the target, B/S and detector are all at the same temperature. 

3.4. 2.2 Calibration Equations 

/ Targ interferogram recorded when viewing a target of interest at temperature = T Talg 

I Co^ interferogram recorded when viewing a cold blackbody calibration target at 
temperature = T CoId 

I w a rm interferogram recorded when viewing a warm blackbody calibration target at 
temperature = T Warm 

U frequency in wavenumbers (cm -1 ) 

c Tar g ( w ) FF7(/ Targ ) = uncalibrated complex spectrum of the target of interest at 

temperature = 7 Targ 

c Coid( v ) FFT(I Cold ) = uncalibrated complex spectrum of the cold blackbody calibration 
target at temperature = T Cold 

c warm( v ') FFT(I Warai ) = uncalibrated complex spectrum of the warm blackbody 
calibration target at temperature = T Warm 

0 instrument phase = the small angle of rotation of the instrument emission spectrum 
from 180 ° 

r(-u) complex instrument spectral response 
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R(u) scalar instrument spectral response 

I Instr interferogram due to instrument emission = the instrument (self-emission) 
contribution to / Targ 

B i ns t r (u) spectral radiance of the instrument (self-emission) in W cm" 2 ster'/cm 1 

^Targ( u ) spectral radiance of the target of interest in W cm 2 ster '/cm"' at temperature 

= T 

■‘Targ 

Bcoid(' u ) spectral radiance of the cold blackbody calibration target in W cm 2 ster'/cm 1 = 
the Planck function at temperature = T Cold 

B warm (ti) spectral radiance of the warm blackbody calibration target in W cm 2 ster '/cm 1 
= the Planck function at temperature = T Warm 

T Instr the temperature of the CIRS instrument - 170.0 K 

The “target” can either be a nominally unknown science target, the “cold” target (deep 
space at 3K), or the “warm” on-board target at the optics temperature, nominally 170K. 
Following (Carlson, 200s) defining c(u) to be the uncalibrated spectrum (FFT of the 
interferogram I, c = FFT(I) ), and proceeding from E^. 16 above, 


Ctarg = r(u) ( (B targ - Bb/s) + (Bdet - B B /s) exp(-iA)). (17) 

A is the constant 71 in the case of a dielectric; in the case of a thin metal film B/S 
(Brasunas, 2002) A has a term linear in wavenumber, indicating a relative shift of the 
interferogram corresponding to the 2 rld term in (17), compared with the first term. 

We then derive calibration equations identical to those of (Carlson. 200s) : 


r(u) = FFT ( I warm - I C oid) / (B warm - B co i d ), and (18) 

B targ = FFT (I targ - I co i d ) / r(u), ( 1 9) 

ignoring B c0 | d (B at 2.7K). These are the calibration equations for FP 3 and 4. 

For FPi, the isothermal instrument means Fb/s = Fdet, and (16) becomes: 


dF mo d = 2 RT(F far - F det ) cos 3 . (20) 

For an external blackbody target at the optics temperature T ins t (nominally 170K for FPi), l ta rg 
is identically zero. Thus the FPi version of Ecp 18 is: 
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I"(u) — FFT ( I cold ) / (Bjnstr B c0 |d) — FFT ( I cold ) / Bjnstr 


(21) 


the same as f Carlson, 2005) . The equation for B tar g is the same as (19). 

3.5 Post processing: resampling and apodization 

After the calibrated power spectrum has been computed, with radiometric and 
wavelength calibration, two further processes are applied that affect (i) the spectral 
resolution and (ii) the sampling grid. In fact, these are applied during the Fourier Transform 
process that creates the power spectrum, to avoid unnecessary repetition of forward and 
reverse transforms, but for simplicity we will present these as independent additional 
processes. 

3.5.1 Apodization) 

Apodization (literally ‘removal of feet’) is a process that smooths the power spectrum, 
dampening down ripples/ringing that are due to the bandwidth limitations of the finite 
Fourier Transform, at the expense of broadening spectral lines. A detailed discussion of this 
topic is beyond the scope of this User Guide (see for example Section 13.4.1 of (Press, 2007)), 
but a short description is necessary, as the user needs to understand how to model the data. 
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Comment [ 1 ]: I don’t think it is necessary to 
show here all possible apodization functions. I 
would recommend discussing only those relative 
to actual CIRS database: boxcar and Hamming. I 
also recommend that the function F(v) ais 
explicitely given for Hamming case (as a function 
of v = frequency (cm' 1 ) and Av = Rayleigh 
parameter . 
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Figure 19: apodization functions in the spatial and spectral domains 

Figure 19 shows the two apodization functions used by CIRS. In the spatial domain an 
apodization functions f is applied to the interferogram by pointwise multiplication. In the 
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spectral domain, this convolves the spectrum with a smoothing function F (columns 2&4). 
So: 

S*F = FFT{Ixf] (22) 

where f and F are a Fourier pair. The boxcar function corresponds to the ‘no apodization’ 
case, and results in a sharply-peaked spectral line (the Instrumental Line Shape, or ILS, is a 
sine function, (sinx)/x with significant side-lobes, or ‘ripples’. However application of the 
Hamming apodization function reduces the side-lobes very substantially, at the cost of 
widening the main lineshape. Quantitatively the FWHM of the Hamming-apodized ILS is 
(0.907/0.603) or 50% greater than for the sine function, while the distance to first ‘null’ (x- 
axis crossing of the ILS) has doubled. This latter quantity is sometimes loosely referred to as 
the Rayleigh definition for spectral resolution, although Rayleigh’s definition technically 
referred to diffraction-limited imaging (when the maximum of the image of one point source 
falls on the first minimum of the second point source). 

In the spatial domain Hamming apodization is applied by the formula: 

f(x) = 0.54 - 0.46 cos (“) (23) 

where a is the maximum path difference. CIRS spectra are computed for unapodized 
(boxcar) and apodized (Hamming) interferograms: however only the apodized spectra are 
provided on the PDS volumes (in DATA/APODSPEC) as apodized spectra are almost 
exclusively used for science. The user must convolve any spectral models accordingly (see 
the ‘Nyquist’, ‘FWHM’ and ‘Rayleigh’ fields for information on spectral resolution). The 
eanvolutite function, which is the F ast Fourier T ransform of the Hamming function, is given 
by: 


F{v) = 


.08 - 0 . 64 a 2 v 2 jsinc( 2 ^:av) 
l-4flV 


(24) 


The respective Rayleigh resolutions are: 

Av=l/ 2 a 

(Boxcar) 

( 25 ) 

And 

Av=l/a 

(Hamming) 

(26) 


Figure 20 shows an example of the effect of apodization on a typical Titan spectrum. 
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Figure 20: example of Titan T23 FIRNADCMP north polar spectrum (average of 800 FP3 spectra), showing the effect of 
apodization (blue) versus the unapodized spectrum (red). 

3.5.2 Resampling 

Whereas the spectral bandpass of each focal plane is set by the interferogram data 
sampling rate (after resampling), the spectrum point spacing is the reciprocal of twice the 
length of the scan (1/2L’), where L’ is measured after ZPD. This typically results in a non- 
integer native point spacing (the ‘Nquist’ value in the database). For ease of modeling 
therefore, the data are resampled to integer values (or simple fraction thereof) grid. To 
avoid loosing information, the resampling grid must be at least as dense as the resolution 
(e.g. i/L’ in the unapodized case). 

As previously mentioned, Fast Fourier Transforms typically use arrays sizes that are 2 n in 
size, therefore simply zero-padding the interferogram to the next larger 2 n array size rarely 
gives the desires point spacing either. Therefore, a two-step process is used. 

1) the interferogram is zero-padded by a large amount (typically a factor - 8 ), to 
create a resulting over-sampled spectrum on an ( 8 x) finer grid than previously, 
then: 
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2) the over-sampled spectrum is interpolated onto the exact wavenumber grid 
desired. Linear interpolation is used, to avoid introducing unwanted artifacts that 
could be caused by use of higher-order functions. 

Table 16 lists the native and resampled point spacings in the CIRS dataset. 


Table 16: native and resampled spectrum point spacings 


RTI (scan length 
in 1/8 s units) 

Native point 
spacing (‘Nquist’) 
(cm 1 ) 

Resampled 
Point Spacing 
(cm 1 ) 

388(39 

7.96 

5.00 

40 

6.53 

5.00 

56 

3-25 

2.50 

96 8( 97 

1.39 

1.00 

224 &. 225 

0.497 

0.250 

400 &. 401 

0.263 

0.250 


3.6 Estimation of Noise Equivalent Spectral Radiance (NESR) 

The sensitivity of the CIRS instrument is specified by the Noise Equivalent Spectral 
Radiance (NESR). The NESR is the target radiance level for which the instrumental signal to 
noise (S/N) ratio is unity. The NESR is defined as: 

NESR(d) = NEP(u)/h Sys Ti 0pt iloiff AQ Ad SQRT(t Eff )] (27) 


where 

- NESR(d) is the Noise Equivalent Spectral Radiance in W cm’ 2 sr 1 / cm 1 


- NEP(d) 

is the noise equivalent power of the detector in W Hz 

-d 

is the frequency in wavenumbers (cm 1 ) 

- *lSys 

is the system efficiency 

" 'lOpt 

is the optical efficiency 

-IrClDiff 

is the loss due to diffraction 

-A 

is the telescope area (cm 2 ) 

-Q 

is the solid angle of the field of view (sr) 

- Ad 

is the spectral resolution (cm 1 ) 

- t E ff 

is the effective integration time (sec) 
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The instrument has the capability for any unapodized spectral resolution from 0.25 to 
10.0 cm" 1 , by commanding the scan time. The apodized spectral resolutions 0.53, 1.00, 2.85, 
6.60, and 15.67 cm 1 have been selected as representative base resolutions for carrying out 
the majority of the CIRS science observations. The measurement parameters for these CIRS 
observing modes and their associated integration times are given in Table 17. 


Table 17: spectral resolution, scan time (RTI), effective integration and frame times 


Unapodized 

Apodized 

Required 

^Frame 

tEff 

tDead — 

Resolution 

Resolution 

RTI 

(sec) 

(sec) 1 

t Frame “ tEff 

= 1/4 U (cm 1 ) 3 

= 1/2 L 2 (cm 1 ) 

(Vs sec) 



(sec) 

0.27 

0.53 

400 

50.000 

47.10 

2.90 

0.50 

1.00 

224 

28.000 

26.11 

1.89 

1.43 

2.85 

96 

12.000 

11.00 

1.00 

3.30 

6.60 

56 

7.000 

5-93 

1.07 

7 .8 5 2 

15-67 2 

38 

4 - 75 ° 

3.10 

1.65 


! t Eff = Effective integration time 
Symmetrical two-sided interferogram 

3 L 2 = scan length from ZPD to the end of scan (cm) ; resolution is here defined as distance from peak to first null 
of the instrument function (Rayleigh definition) 

The calculation of the effective integration time (tEff) is discussed in Appendix D. The frame 
time is the repeat cycle for the scan and is approximate. 


3.6.1 NESR Conversions 
Spectral Resolution 

The NESR for a given spectral resolution (NESR 2 ) may be calculated from the NESR for the 
reference resolution (NESR, = 0.53 cm 1 ) by: 

NESR, = (Av,/Av 2 )*[SQRT(tEf f1 /t Ef f 2 )]*NESR, (28) 

Using Table 17, the conversion factors are: 

NESR(i.oo cm 1 ) = 0.7118 * NESR(o. 53 cm' 1 ) (29a) 
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NESR(2.85 cm' 1 ) = 0.3848 * NESR(o.53 cm' 1 ) 


(29b) 


NESR(6.6o cm 1 ) = 0.2263 * NESR(o.53 cm' 1 ) 


(29c) 


NESR(i5-67 cm' 1 ) = 0.1318 * NESR(o.53 cm' 1 ) 


(29d) 


Number of Pixels 

For FP3 and FP4 scans, individual pixels may be summed together. The improved NESR 
would then be: 


For example, if the detectors are used in PAIR mode then the NESR is the single pixel 
NESR/n/2. 

3.6.2 NESR Spreadsheets and Plots 

The CIRS NESRs were calculated from 100 scan averages of deep-space interferograms 
recorded during the Jupiter encounter on days 330 and 332, 2000. Each NESR is apodized 
and corresponds to a single scan. The measurement conditions were: 

Target temperature = 2.7 K (deep space) 

FPi: Far-IR polarizer temperature = 170.04 K 

FP3 8( FP4: average optics assembly temperature = 169.84 K; FPA temperature = 77.74 K 

NESRs for apodized spectral resolutions 0.53, 1.00, 2.85, 6.60, and 15.67 cm' 1 are 
tabulated for focal planes 1, 3, and 4 in files provided in the CALIB directory of the CIRS PDS 
archive volumes (CIRS_NESR_FPn.ASC, where n=i,3,4). Plots of the FPi, FP3, and FP4 NESRs 
for each spectral resolution are shown in Figs. 21-24 . 


NESR(N pixels) = [i/SQRT(N)]*NESR(i pixel) 


(30) 
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Figure 21: NESR for FPi versus spectral resolution 
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Figure 22: NESR for FP3 versus spectral resolution 
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0.53 cm ' 1 Frame Time = 50.0 sec. 

1 .00 cm ' 1 Frame Time = 28.0 sec. 

2.85 cm ’ 1 Frame Time = 12.0 sec. 

6.60 cm ’ 1 Frame Time = 7.0 sec. 

15.67 cm ' 1 Frame Time = 4.75 sec. 
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Figure 23 : NESR for FP 4 versus spectral resolution 

3.7 Discussion: limitations of standard spectra products 


Several final notes and caveats are given here about the limitations of the standard 
spectral data products, above and beyond previous discussions of noise, interferences etc. 

3.7.1 Bandpass and aliasing 

CIRS is a Fourier Transform Spectrometer (FTS), so the bandpass of each focal plane is 
set by the sample spacing of the native interferogram. The highest wavenumber uniquely 
measured is vNY=(o.s/AL) cm' 1 , where AL is the sample spacing in cm. For example FPi 
acquires -2800 samples over -1 cm of mirror travel (2 cm of light path difference), so 
AL=o. 000714 cm and the Nyquist frequency (wavenumber) is -700 cm 1 . Light at higher 
wavenumbers v’ is aliased into the FPi band at frequency v, where: v=2v NY -v’. E.g. 1000 cm 1 
is the alias of (2x700 - 1000) = 400 cm 1 . To avoid the uncertainty caused by aliasing, a filter 
must be applied to block the unwanted higher frequencies: CIRS achieves this in the receiver 
electronics by numerical filtering. 

FP3 and FP4 are similar to FPi, except that they are designed to operate in the second 
(v ny - 2 v ny ) and third (2v NY -3v NY ) orders respectively, so must be filtered at both the high and 
low ends of the bandpass. In practice, no filter F(v) has an infinitely-sharp cut-off, rather 
there is a gradual roll-off. .Figure 24. shows the high-frequency end of the FP4 range. This 
gradual cut-off requires us to have some criterion for defining the acceptable ‘clean’ spectral 
range. The CIRS definition is such that F(v’)/F(v)<i% - i.e. the spectral range where the ratio 
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of the weighting of the aliased frequency (v’) is less than 1% of the in-band frequency (v) 
weighting. 



Table 18: actual and clean bandpasses for each focal 
plane 


FP 

Signal 

Clean 


bandpass 

Bandpass 

1 

0-708 

0-651 

3 

556-1133 

604-1094 

4 

1020-1529 

1131-1433 


1300 1350 1400 1450 1500 1550 1600 1650 1700 

wavenumber (1/cm) 

Figure 24: frequency filter applied to short wavelength end of 
CIRS FP4. 


The nominal and clean bandpasses for each focal plane are given in Table 18. One point of 
particular concern is that, while the ‘clean’ ranges of FPi and FP3 overlap, those of FP3 and 
FP4 do not: meaning that there is a spectral range from 1094-1131 cm 1 where aliasing may 
become significant. Fiowever, any signal seen at, for example, 1110 cm 1 on FP3 may be 
checked for reality by looking at the same frequency on FP4. More problematic is the range 
from 1433-1529 cm" 1 on FP4, as there is no second measurement with which to compare the 
spectrum. The only solution is to model the aliasing: to extend the range of the spectral 
calculation to the highest aliased frequency, and then to combine the nominal and aliased 
ranges with the appropriate scaling F(v) : 

i’(v) = [F(v)l(v)]+[F(2vNY-v)l(2vNY-v)] (31) 


Figure 25 shows an example Titan spectrum, with the ‘clean’ bandpasses shown in blue and 
the compromised bandpass in red. The aliased bandpass of FP4 includes the vs band of C 2 Fl6 
(1468 cm 1 ) and the v 7 (1464 cm 1 ) and v 24 (1472 cm 1 ) emission bands of CjHs. It is an area of 
on-going work at this time to add a model of the aliasing behaviour into spectral synthesis 
models. 


66 


Brightness Temperature (K) 



Figure 25: example CIRS high resolution (0.5 cm' 1 ) spectrum, with 'clean' bandpasses for each focal plane in blue and 
'aliased' ranges in red. 
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4-0 How to search and extract data 


4.1 Data Format Description 

The CIRS PDS dataset is divided into volumes, each of which is a self-contained unit 
comprised of CIRS and ancillary data for one month of the Cassini mission. The volumes are 
composed of data, formally known as ‘Time Sequential Data Records’ (TSDR), related search 
software, and other documentation. A schematic diagram of the creation of the Cl RS archive 
volumes is shown below (Figure 26). 



Figure 26: schematic of CIRS archive volume creation 


The organization of the database is now described in more detail. 
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Figure 27: relative locations of CIRS detectors, showing numbering of pixels and q-points. 


4.1.1 Data Levels 

The level -0 tables (base level) are contained in /DATA/UNCALIBR/. These are described in 
detail in the document DATASIS.PDF. These tables types (IFGM, OBS, I HSK, FRV) are 
essentially the raw-level information which came from the instrument: the only process 
applied was unpacking (uncompressing) the instrument packets and, in the case of the 
housekeeping data, interpolating the 64-sec interval housekeeping data onto the science 
scan intervals. Original housekeeping data is stored in /DATA/FISK_DATA/, which is read out 
at 64-sec intervals when the instrument is either ON or in SLEEP mode. Pointing and 
geometry information is stored in /DATA/NAV_DATA/. These tables are described in detail in 
the DATASIS document, and also in the FMT format files in the NAV_DATA/ area. Primarily, 
these store the spacecraft attitude and position with respect to major bodies (GEO file 
type), the position of the 11 detectors (99 q-points, or fiducial reference marks) on all bodies 
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in the FOV (POI file type, see Figure 27), and the list of bodies seen in each detector (TAR 
type). 

Pointing and geometry information is derived from the spacecraft re-constructed 
attitude ('C') and trajectory ('SPK') kernels, using NAIF (Navigation and Ancillary Information 
Facility) toolkit routines supplied by JPL. Calibrated spectra are stored in two varieties: 
apodized fwith Hamming function) and unapodized. These are stored in /DATA/APODSPEC/ 
and /DATA/U NAPS PEC/ respectively. Apodization is the process of mathematical filtering or 
smoothing, which removes 'ringing' effects from finite-width mathematical FFTs, widening 
the instrument line shape in the process. Either form is valid for science, depending on user 
requirements, although apodized is often preferred as it can be easier to model 

4.1.2 Vanilla database 

The data is in the form of binary tables, also known as the Vanilla database. Vanilla is a 
simple database access tool, also included on each volume as both source and executable. 
Note however that the Vanilla software is NOT REQUIRED in order to read the binary data, 
but it may simplify the task. Vanilla links together all binary files into a large, 'logical' table, 
which can be searched based on key fields. However, the user may read the binary data 
directly, using any programming language desired. If this approach is taken it is extremely 
important to pay attention to binary field widths: field widths vary, and typically trying to 
read an entire record into a data structure will fail: if it assumes padding of fields. The user 
must read the data in field-by-field, specifying width. Field widths are described in the *.FMT 
files, which are placed in each data directory. Some tables have both fixed-length and 
variable-length fields. In this case, the variable length part of the record is stored separately, 
in a file with the same name but an extension .VAR, instead of .DAT. The third type of file 
extension is .LBL, which is the PDS detached label. 

Important note : by convention hereafter we will refer to tables by uppercase names 
(IFGM, TAR, RIN etc) and fields within the tables by lowercase names (scet, det etc). This is 
to avoid the confusion that may arise when the same term is used for both a table type, and 
a field within that table (ifgm, fifm, ispm). 

4.2 Vanilla Search Tool: Overview 

The Vanilla tables were designed to be a relational database, with different tables linked 
by keys (i.e., fields common among some or all of the various tables). There are several 
important advantages to this approach. Being able to store the data in a binary format 
without having to repeat all fields for every possible combination of scet, detector, and 
Saturnian body produces enormous savings in storage space. By storing a given table in 
multiple file fragments, search times are minimized by specifying only the applicable 
fragment or fragments on the vanilla command line. Finally, because the source code is 
available, modifications were possible that allowed customization of the database for CIRS- 
specific needs and for PDS standards compliance. 


Table 19: table keys at a glance 

TABLE NAME KEY 1 KEY 2 KEY 3 
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DIAG 

scet 

det 

FRV 

scet 


GEO 

scet 

body id 

HSK 

scet 

det 

IFGM 

scet 

det 

IHSK 

scet 


ISPM 

scet 

det 

OBS 

scet 


POI 

scet 

det 

RIN 

scet 

det 

TAR 

scet 

det 


An Excel spreadsheet with all fields and brief descriptions is in Vanilla Field List.xls. 


Vanilla is a command line program. This enables it to be easily used within other 
programs and scripts by "piping" its output to the next process. To use it you'll need access 
on your computer to a terminal (or command) window from which you can enter the 
commands, the vanilla executable file, and the vanilla executable command within your 
shell's command search path. 

The vanilla executable file is available on PDS volumes in the SOFTWARE/BIN directory 
and has been precompiled for the Linux and Solaris operating systems. The source code is 
available in the software/src/vanilla directory should you have an OS different from the 
above and want to compile it yourself. Only a single executable file is needed to run the 
vanilla database query program. 

4.2.1 vanilla help 

Typing only the vanilla command name at the command prompt will give you a quick 
syntactical description of how to use the program. It will look similar to: 


Vanilla version 3.4.21 

Usage: vanilla [options] dataset 

dataset: (required) 

directory-name containing the DATASET file -OR- 
dataset-file-name -OR- 

-files data-filei [data-file2 data-file3 ...] 

query options: 

-fields 'fieldi field2 ...' (required) 

-select 'fieldi low high field2 low high ...' (optional) 

-Bfm filename binary output, cooked 4 byte words 
-Bdm filename binary output, cooked 8 byte words 
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index generation options: 
-index 'fieldi field2 


(required) 


index based count options: 

-select 'fieldi low high field2 low high (required) 

-count 

-OR- 

-select 'fieldi low high field2 low high ...' (required) 

-idxstats 

Note: 

Use one type of option at one time, i.e. don't use query with index etc. 


Of these many possibilities, by far the most commonly used are just the "dataset" 
specification and the "-fields" and "-select" query options which means your query line will 
look more like: 

vanilla -fields "fieldi field2 ..." -select "selecti selectimin selectimax ..." dataset_path 
where: 

fieldi... are the fields you would like extracted from the database 
selectX is a field which you would like to restrict to a range 
selectXmin is the minimum value allowed for selectX 
selectXmax is the maximum value allowed for selectX 
dataset_path is the path to the DATASET file 


4.2.2 Database structure 

Each CIRS PDS volume contains extensive documentation on its structure. For a detailed 
description of the CIRS database including its file structure in the volume, the table layout, 
and the data fields available within each table, see the DATASIS.PDF document 

4.2.3 Dataset usage 

The dataset_path mentioned above in the vanilla command description can be specified 
several different ways, and, depending on which you choose, can dramatically shorten the 
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time necessary to complete a query. The dataset_path can either point to a directory with a 
dataset file in it or it can point to a dataset file outside of the PDS volume. A dataset file is 
used to tell vanilla what tables should be included in the search and/or where to search for 
another dataset file. 

For example, the dataset file in the PDS volume's data/apodspec/ directory is 
ISPM 

,./NAV_DATA/ 

This tells vanilla to look in the current directory for the ISPM tables and to look in the 
NAV_DATA directory for another dataset file. The dataset file in NAV_DATA contains 

../UNCALIBR/ 

TAR 

RIN 

POI 

GEO 

and the dataset in UNCALIBR is 

OBS 

IHSK 

IFGM 

FRV 

DIAG 

Therefore, by pointing to the dataset file in the apodspec/ directory, you will be 
searching across all available tables in the database. Flowever, if you only wanted to query 
for fields contained in the OBS table, you could just point to the dataset file in the UNCALIBR 
directory, allowing vanilla to skip the overhead of including the navigation and spectra 
tables. 

If your output fields are restricted to only a single table, there is one more way to greatly 
decrease the number of files to be searched. By using the "-files" option (see the full vanilla 
help output) rather than a dataset specification, you can specify exactly which files to 
search. For example, if you knew your data was completely self-contained within the OBS 
table for 4 January 2008, in the command line instead of dataset_file you would specify 

-files path/OBSo8oi0400.LBL path/OBSo8oi0404.LBL path/OBSo8oi04o8.LBL 
path/OBSo8oi04i2.LBL path/OBSo8oi04i6.LBL path/OBSo8oi0420.LBL 
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where path is the path to the label files of interest. While it is an extremely restrictive case, it 
will result in a very fast search. Another advantage of using this approach is that wildcards 
are allowed. So the above could have been shortened to 

-files path/OBSo8oi04?? 

At the other extreme, suppose you wanted to query the entire database consisting of 
multiple PDS volumes. Assuming you've copied all the volumes to your hard drive at /data, 
then the dataset file would look like 

/data/cocirs_040i/data/apod_spec/ 

/data/cocirs_0402/ data/apod_spec/ 

/data/cocirs_0403/ data/apod_spec/ 


/data/cocirs_i705/ data/apod_spec/ 

/data/cocirs_i7o6/data/apod_spec/ 

Note, when specifying a directory within a dataset file, always place a "/" at the end of the 
path. If you don't, vanilla will output a cryptic "No tables in dataset." message. 


4.2.4 ‘-fields’ parameter 

The fields listed with the fields parameter specify the columns that will be output to the 
screen. Any field listed in the .FMT file for each table can be used. In addition, if that field is 
an array, you can specify a single element using square brackets or a if a range of 
elements is used. For example, ispm[i:io] would extract the first 10 elements of the 
spectrum while spm[gg] would extract only the ggth (vanilla array indices start at 1). 
Typically the entire array can be extracted with just its name, but in the case of a variable 
length array (like ispm), the entire array can only be output if "ispm[]" is used. It does not 
matter whether you use upper- or lower-case to specify the field. The field can be specified 
by either using its field name alone or using both the table name and the field name 
together, separated by a e.g., "OBS.scet". 


4.2.5 ‘ select’ parameter 

The selections listed with the select parameter will restrict the output field columns to only 
those records matching the selection range. Any field listed in the .FMT files for each table 
can be used. In addition, if that field is a fixed-length array, you can specify an element of the 
array, for example, "POI.Iatitude_zpd[5] -10 10". Multiple selection parameters are logically 
"AN D'ed" together, however, if two or more of the selection fields are the same (e.g., when 
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you want multiple, disjoint scet ranges) those fields are logically "OR'd". Like the fields 
parameter, case does not matter, and the field can be specified either alone or with its table 
name. 


4.2.6 vanilla output 

By default, all vanilla output is ASCII text written to screen with tab characters separating 
output fields and a return character at the end of each line (i.e. a "record"). The very first 
line of output will always be a listing of the field names requested. 


4.2.7 Tips and tricks 

These are tips that may help when creating vanilla queries. Some have already been 
mentioned above but are repeated here because of their usefulness. 

1. A directory path in a dataset file ALWAYS needs a trailing "/". 

2. Double quotes (" ") or single quotes (' ') can be used to enclose the fields and/or the 
selections on the vanilla command line. 

3. The dataset_path, fields, and select parameters can appear in any order on the 
command line. 

4. You CANNOT specify a single element from an array in the selection parameter of the 
vanilla command line if that array has variable length. In short, this means you cannot 
use single elements from the calibrated spectrum (ispm) or interferogram (ifgm, 
fifm) arrays in the select parameter. For example: 

vanilla myDatasetPath -fields "fieldi field2..." -select "ispm[gg] 5e-8 6e-8" 
will not work! 

5. To shorten your query time, if you need data from only a single month or a subset of 
months, create a dataset file with only those months listed. 

6. If the field or selection you want is an entire array, you can specify it just with its 
name. But if it's a variable length array, square brackets must be included. The only 
variable length field arrays are currently frv, ifgm, and ispm. 

7. To prevent confusion, it is recommended that a field that is also a key be qualified 
with its table name, e.g., OBS.scet. 
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8. Note that scet, as used in the tables, typically implies when the measurement was 
taken, at what time the value was calculated, or with what time the field should be 
associated (usually that of an IFM). For example, IFGM.scet indicates the time at the 
start of the interferometer's scan, POl.scet gives the time at which the pointing 
information is applicable (specifically chosen to be that of an interferogram), and the 
scet time for DIAG.tzpd indicates which IFM it applies to. 

9. The tool used to query the CIRS data set was designed to be a generic, simple, lean 
program for querying data with few frills. Thus the name "vanilla" was chosen. 

10. Vanilla is case insensitive with the “fields” and “select” query strings (enclosed in “”), 
although the directory/dataset file name is case sensitive on UNIX operating systems. 


4.3 Vanilla and the CIRS Database 

Another of the great advantages of a relational database is that each table does not 
need to be one to one with the other tables, i.e., they may not all have the same number of 
records. Flowever, to be useful, a large number of their keys should match across tables. It is 
through the concept of "keys" that a relational database can match records from one table 
to those in another. Understanding how keys are implemented within the CIRS database is 
crucial to creating meaningful queries with the intended output. 

The primary key field, used to cross-link CIRS data in different files together, is spacecraft 
event time (scet) stored as an integer number of seconds since 1/1/1970 00:00:00 UT, 
ignoring leap seconds. Also stored as retrievable fields, but not key fields, are spacecraft 
clock (sclk), and the sclk partition number. These are stored as two integers. 

sclk is essentially an integer number related to oscillations of the on-board quartz crystal 
clock. Partition number is advanced by one when the sclk count is reset, expected to be a 
rare event, if ever, sclk ticks are not even linearly related to universal time, and the 
conversion between the two requires a polynomial to account for drifts in the sclk oscillation 
frequency, scet on the other hand is linearly related to universal time. 

scet is used as the primary time key field rather than sclk for two reasons: firstly, sclk 
may reset, and not be monotonically increasing, in which case the partition number also 
changes; and secondly, scet may easily and conveniently be converted to text time using 
standard computer library functions, whereas interpreting sclk requires the NAIF toolkit and 
Cassini time kernel (SCLKSCET) files. 

At the time of writing, standard computer library routines are available which facilitate 
the conversion to and from text time to scet number, such as the C and Perl functions 
gmtime and timegm. 
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The golden rule that vanilla will use when searching through records from different 
tables is simple - when the selection criteria are satisfied, IT WILL ONLY COMBINE RECORDS 
IF THE COMMON KEYS AMONG THE TABLES CAN BE MATCHED. A key or keys within a table 
are used to uniquely identify a record (consisting of those keys plus all the other fields in the 
table). For tables using only a single key, each key value is unique. For tables with multiple 
keys, their key combinations are unique. Some CIRS-specific examples should help make this 
clear. (Note that because of improvements we continually make to our database, you may 
not be able to exactly replicate these results.) 

Example i: 

Query a single table (OBS) for a single field given a selection criteria on scet. 

vanilla /data/cocirs_050i/data/apodspec/ -fields "OBS. scet" -select "OBS. scet 1104538500 
1104538542" 

yields 

obs.scet 

1104538500 

1104538506 

1104538542 

The OBS table has only 1 key ("scet"). This output implies there are 3 records in the OBS 
table that match the selection criteria. Of course, if any of the other OBS fields had been 
included in the query, they would have been output as well. 

Example 2: 

Query a single table (ISPM) for two fields given a selection criteria on scet. 

vanilla /data/cocirs_050i/data/apodspec/ -fields "ISPM. scet ISPM.det" -select "ISPM. scet 
1104538500 1104538542" 

yields 

ispm.scet ispm.det 
1104538500 0 

1104538500 11 

1104538500 13 

1104538500 15 

1104538500 17 

1104538500 19 

1104538500 22 

1104538500 24 

1104538500 26 
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11045385°° 28 
11045385°° 3 ° 

1104538542 o 

Here it is clearly evident why ISPM needs two keys to keep each key combination unique. 
If the table had been designed with only scet used as a key, multiple records for the same 
time but different detector numbers would have the same key (“degeneracy”). 

Example 3: 

Query fields from two tables given a selection criteria on scet. 

vanilla /data/cocirs_050i/data/apodspec/ -fields 'OBS.scet OBS.rti ISPM. scet ISPM.det 
ISPM.ispts' -select "ISPM. scet 1104538500 1104538542" 

yields 


obs.scet 

obs.rti ispm.scet 

ispm.det 

1104538500 

38 

1104538500 

0 

139 

1104538500 

38 

1104538500 

11 

112 

1104538500 

38 

1104538500 

13 

112 

1104538500 

38 

1104538500 

15 

112 

1104538500 

38 

1104538500 

17 

112 

1104538500 

38 

1104538500 

19 

112 

1104538500 

38 

1104538500 

22 

95 

1104538500 

38 

1104538500 

24 

95 

1104538500 

38 

1104538500 

26 

95 

1104538500 

38 

1104538500 

28 

95 

1104538500 

38 

1104538500 

3 ° 

95 

1104538542 

38 

1104538542 

0 

139 


Comparing the output of examples 1 and 2 shows that the time for scet = 1104538506, 
from Example 1, is not listed. It has been excluded because there is no matching scet in the 
ISPM table. This is an example of the golden rule in practice - records from both tables are 
combined only if the values of their common keys match (the only common key here is scet). 

Example 4: 

Query a single table with a key field, THE TABLE UNDESIGNATED, and a scet selection, 
vanilla /data/cocirs_050i/data/apodspec/ -fields 'scet' -select "scet 1104538500 1104538542" 
yields 
scet 

1104538500 
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11045385°° 

11045385°° 

11045385°° 

11045385°° 

11045385°° 

11045385°° 

11045385°° 

11045385°° 

11045385°° 

11045385°° 

1104538542 

This example illustrates how vanilla handles ambiguities that arise when the requested 
field occurs in more than one table. The query above returned the scets from the ISPM 
table. The duplicate scets are returned because using scet alone is not sufficient to identify a 
unique record; in the ISPM table, both scet and det (the table key set) are required. If vanilla 
encounters a field common to more than one table, it will search the list of tables and 
directories listed in the dataset file and use the first table that has that field. In this case, 
ISPM is the first table in the dataset list and therefore ISPM. scet is used. If you want to avoid 
this type of ambiguity and prevent possible confusion in your query's output, just USE THE 
TABLE NAME WITH THE FIELD WHEN THAT FIELD OCCURS IN MULTIPLE TABLES. In most 
practical cases, this means the fields "scet" and "det" should be qualified with the table 
name in your queries. 

4.4 Vanilla extraction of pointing/geometry 
NAVIGATION TABLE QUERIES, GEO: 

Of the navigational tables, the RIN and TAR tables both use scet and det fields as the 
keys to identify a unique table record. However, the GEO and POI tables each use a very 
different key set and require special consideration when used in vanilla queries. 

The GEO table was designed to provide information on the positions and orientations of 
the local object bodies and the spacecraft. As seen in Table 19, the field keys for the GEO 
table are scet and bodyjd. The field "det" is not needed as a key, nor even as a field in the 
table, because the relative positions between the spacecraft and a body are independent of 
detector. The bodyjd field is a subset of the bodies in the Jupiter and Saturn systems with 
values defined by the NAIF ID for those bodies (see Table 3 "Selected NAIF IDS of relevance 
to Cassini CIRS" in the DATASIS.PDF PDS document). The GEO table's primaryjd field 
identifies which system is relevant for that scet. For each scet, information for all bodies 
within the primaryjd field's designation will be available. Therefore, for the relatively brief 
period that Jupiter was encountered, each scet will have entries for body id's of 501-516 and 
599. Likewise, for the prime mission and beyond, the primaryjd field will be 699 (Saturn) 
and each scet will have entries for bodyj'd's of 601-618 and 699. This means at Saturn, each 
scet in the GEO table will have 19 entries! Not specifying a bodyjd when using fields from 
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the GEO table is a common cause of confusion as many more records will be produced than 
expected with, seemingly, duplicate scet times. Example 5 demonstrates this. 

Example 5: 

Query for the GEO key set fields using a selection criteria on scet. 

vanilla /data/cocirs_050i/data/apodspec/ -fields 'GEO. scet body_id' -select "GEO. scet 
1104538500 1104538506" 

yields 

geo.scet bodyjd 
1104538500 601 

1104538500 602 

1104538500 603 

1104538500 604 

1104538500 605 

1104538500 606 

1104538500 607 

1104538500 608 

1104538500 609 

1104538500 610 

1104538500 611 

1104538500 612 

1104538500 613 

1104538500 614 

1104538500 615 

1104538500 616 

1104538500 617 

1104538500 618 

1104538500 699 

1104538506 601 

1104538506 602 

1104538506 603 

1104538506 604 

1104538506 605 

1104538506 606 

1104538506 607 

1104538506 608 

1104538506 609 

1104538506 610 

1104538506 611 

1104538506 612 

1104538506 613 
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110453 8 5°6 614 

1104538506 615 

1104538506 616 

1104538506 617 

1104538506 618 

Another important concept within the GEO table is that of timing. In order for different 
tables to be joined productively, the majority of the scets within each table should be the 
same. The scet for the ISPM table refers to the first recorded raw sample of the 
interferogram. Values of fields in other tables for scet are then either for measurements 
taken at that scet time or interpolated to it (as in the IHSK table). However, in the GEO table 
it is often more important to know the location of an object during the period at which the 
interferometer is at Zero Path Difference (ZPD) because during this time the great majority 
of the signal is being collected. So while GEO. scet indicates to which interferogram the GEO 
field should be associated, the actual positions are calculated at the time of ZPD, typically a 
few seconds later than the scet. In this way the functionality of the database is maintained 
while providing the highest degree of navigational accuracy. 

NAVIGATION TABLE QUERIES, POI: 

The POI tables store the pointing for each active detector at the ZPD of the IFM. One POI 
record is produced for each active detector and each non-rings non-primary target body in 
the field of view plus a record for the primary body (either Saturn or Jupiter) whether or not 
it is in the field of view. The POI table therefore requires 3 fields within its key set to uniquely 
identify a record: scet, det, and targetjd. There are typically 11 channels (detectors) active 
during a given time. The number of targets within the field of view varies from observation 
to observation and therefore so does the number of records with a given scet. Example 6 
displays the available, unique records for the scet range used in Example 5. 

Example 6: 

vanilla /data/cocirs_050i/data/apodspec/ -fields 'POI. scet POI. det POl.targetjd' -select 
"POI. scet 1104538500 1104538506" 


yields 

poi.scet poi.det poi. targetjd 


1104538500 

0 

608 

1104538500 

0 

699 

1104538500 

11 

608 

1104538500 

11 

699 

1104538500 

13 

608 

1104538500 

13 

699 

1104538500 

15 

608 

1104538500 

15 

699 

1104538500 

17 

608 


81 



1104538500 

17 

699 

1104538500 

19 

608 

1104538500 

19 

699 

1104538500 

22 

608 

1104538500 

22 

699 

1104538500 

24 

608 

1104538500 

24 

699 

1104538500 

26 

608 

1104538500 

26 

699 

1104538500 

28 

608 

1104538500 

28 

699 

1104538500 

30 

608 

1104538500 

30 

699 

1104538506 

0 

608 

1104538506 

0 

699 

1104538506 

12 

608 

1104538506 

12 

699 

1104538506 

14 

608 

1104538506 

14 

699 

1104538506 

16 

608 

1104538506 

16 

699 

1104538506 

18 

608 

1104538506 

18 

699 

1104538506 

20 

608 

1104538506 

20 

699 

1104538506 

22 

608 

1104538506 

22 

699 

1104538506 

24 

608 

1104538506 

24 

699 

1104538506 

26 

608 

1104538506 

26 

699 

1104538506 

28 

608 

1104538506 

28 

699 

1104538506 

30 

608 

1104538506 

3° 

699 


The output shows that there are 11 different vanilla detectors per scet, as expected, plus 
an entry for the primary target body (Saturn) which may or may not be in the field of view, 
and for lapetus (targetj'd = 6o8=NAIF ID) which is in the field of view. 

NAVIGATION TABLE QUERIES, TAR: 

To find out if Saturn is in the field of view, you can use the TAR table fields, TAR. saturn 
and TAR.iapetus. They will be one if in the field of view and zero if not (Example 7). 
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Example 7: 

vanilla /data/cocirs_050i/data/apodspec/ -fields POl.scet POl.det POl.targetjd saturn 
iapetus' -select "POl.scet 1104538500 1104538506" 


yields 

poi.scet poi.det poi.targetj'd saturn iapetus 


1104538500 

0 

608 

0 

1 

1104538500 

0 

699 

0 

1 

1104538500 

11 

608 

0 

1 

1104538500 

11 

699 

0 

1 

1104538500 

13 

608 

0 

1 

1104538500 

13 

699 

0 

1 

1104538500 

15 

608 

0 

1 

1104538500 

15 

699 

0 

1 

1104538500 

17 

608 

0 

1 

1104538500 

17 

699 

0 

1 

1104538500 

19 

608 

0 

1 

1104538500 

19 

699 

0 

1 

1104538500 

22 

608 

0 

1 

1104538500 

22 

699 

0 

1 

1104538500 

24 

608 

0 

1 

1104538500 

24 

699 

0 

1 

1104538500 

26 

608 

0 

1 

1104538500 

26 

699 

0 

1 

1104538500 

28 

608 

0 

1 

1104538500 

28 

699 

0 

1 

1104538500 

3° 

608 

0 

1 

1104538500 

3° 

699 

0 

1 

1104538506 

0 

608 

0 

1 

1104538506 

0 

699 

0 

1 

1104538506 

12 

608 

0 

1 

1104538506 

12 

699 

0 

1 

1104538506 

14 

608 

0 

1 

1104538506 

14 

699 

0 

1 

1104538506 

16 

608 

0 

1 

1104538506 

16 

699 

0 

1 

1104538506 

18 

608 

0 

1 

1104538506 

18 

699 

0 

1 

1104538506 

20 

608 

0 

1 

1104538506 

20 

699 

0 

1 

1104538506 

22 

608 

0 

1 

1104538506 

22 

699 

0 

1 

1104538506 

24 

608 

0 

1 
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1104538506 

24 

699 

0 

1 

1104538506 

26 

608 

0 

1 

1104538506 

26 

699 

0 

1 

1104538506 

28 

608 

0 

1 

1104538506 

28 

699 

0 

1 

1104538506 

30 

608 

0 

1 

1104538506 

30 

699 

0 

1 


So Saturn was not in the field of view and was only included because it was the primary 
target body. 

The distinction between the fov_targets field in the TAR tables, and the individual bytes 
for each target can also be useful, since it allows searches that isolate times when two 
objects are in the field of view, or times when only one object is in the field of view. For 
example, the search: 


vanilla . -fields “tar.scet tar.det” -select “titan 1 1” 


will find all entries in which Titan is in the field of view, while 
vanilla . -fields “tar.scet tar.det” -select “fov_targets 4096 4096” 
finds all entries in which only Titan is in the field of view. 


4.5 Final Words On Vanilla 

Additional material on the vanilla database search engine and its use with the CIRS 
database are available within each CIRS PDS volume. The file vanilla-guide.pdf in the 
directory software/doc/ gives a description of the vanilla search engine as well as a detailed 
explanation of the generic database file format required by vanilla. In the same directory, 
the file cirs_vanilla_examples.pdf expands on the examples above, providing further 
guidance on querying the CIRS database using more detailed queries. 

4.6 Other Extraction Tools 

The CIRS team uses a variety of scripts and programs that either supplement or directly 
operate on vanilla output. Their functions include conversions between scet (i.e., UNIX time) 
and a human-readable UTC time, plus a utility to average spectra returned from a vanilla 
query. While currently only being utilized within the team, based on user interest, they may 
be provided on PDS volumes at a later date. 
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5-0 End-to-end example of data use 


In previous sections of this document, we have described in detail the background about 
CIRS and diverse topics about the instrument, observation designs, operations, calibration 
and how to extract data from the binary records. In this section we finally come to the core 
purpose of this CIRS User’s Guide, a detailed discussion of how you - the reader and 
hopefully scientific user - can select the data you want, subset the data in different ways, 
comprehend the limits to the data and instrument modeling specifics, and finally reverse 
engineer the infrared spectrum to find scientific data regarding the target of interest. This 
section begins with an overview of the tasks required and their ordering, to achieve this 
goal. 

5.1 Overview and general algorithm 


A high-level overview of the steps involved from the very beginning of contemplating what 
CIRS data might be useful for a particular scientific investigation, to acquiring the final 
results of the study, are as follows: 

1) Selection of observations: beginning with target, determine types of observations 
available, and then research specific instances of those observations to find the most 
appropriate. 

2) Subset spectra: beginning with all potentially relevant spectra, list their 
characteristics including geometrical parameters, and exclude or subset the spectra 
as required. 

3) Bad data check: check the subseted spectra more rigorously for bad data. 

4) Average spectra: similar spectra are added together, if required, to increase signal-to- 
noise ratio (SNR), priorto modeling. A detailed assessment of SNR forthe average is 
made at this time. 

5) Determine field-of-view model: if the finite field-of-view of the detector convolves an 
inhomogeneous area, then the user needs to consider whether to explicitly model 
this (e.g. by inclusion of multiple ray paths), or estimate it as a systematic error. 

6) Spectral modeling and inversion: the final stage is to model the spectrum, adjusting 
parameters iteratively until a good fit is made. At this point the final parameters 
values are set and error bars estimated. 

The following sub-sections will describe each of these topics in more detail. 

5.1.1 Selection of observations 

In Section “2.2 Science Objectives” the broad scientific objectives of CIRS are 
summarized, while in Section “2.3 Observation Strategies” the actual implementation of 
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these objectives in terms of planned observations are described in depth. If you are new to 
using CIRS data, you should familiarize yourself at this point with the various observation 
types, and determine which you are interested in studying. 

Next comes the issue of finding exactly when these observations were conducted, and 
what the pointing and other parameters were. All CIRS observations implemented and 
archived are listed in a Time Ordered Listing (TO L) file, which can be found in the INDEX 
directory on the PDS archive volumes. The TOL files have a name such a C23_TOL.ASC or 
S49_TOL.ASC, where the first three characters of the name correspond to the spacecraft 
execution sequence, simply a block of time for which the spacecraft executes a pre-set 
series of maneuvers. The earliest sequences (prior to SOI) had names beginning with ‘C’ for 
cruise, while those after SOI had names beginning with ‘S’ for Saturn. The TOL files were 
created by direct export from the Cassini Information Management System (CIMS), an on- 
line database available during the Cassini Mission to all teams, and keep track of all planned 
and executed observations and maneuvers by all instruments and engineering teams. 

The TOLs are text files of a specific text-table format, known as CSV, for comma 
separated values. Each line in the file, known as a record, consists of an identical number of 
data fields, delimited by double quotes, and separated from the other fields by a comma. 
Together all the records form a table, which can be read by standard spreadsheet programs 
such as Excel and Open Office Calc. The TOL records contain the start time and end time of 
each observation in a format yyyy-doyThh:mm:ss, where yyyy=year, doy=day of year 
(beginning with Jan 1 st =001), hh=hour, mm=minute, ss=second. The observation name is 
also recorded. 

Therefore, either by reading the TOLs into a text editor and using standard search 
commands, or by reading into a spreadsheet program, the entire observational database 
may be searched for specific instances of any observation type, using the 9-character stem 
of the observation instance, such as FIRNADMAP or MIDIRTMAP. 

The next task is to convert the DOY-style date and time into a SCET time, as described in 
Section “4.0 How to search and extract data”. Finally, the Vanilla software tool may be used 
to query the database for records between the start and end scet time of the observation, 
to extract pointing information such as latitude, longitude, local time, altitude at FOV center, 
and other parameters. 

5.1.2 Subset observations 


The user should begin by initially extracting all the spectra for the desired observation(s), 
including necessary pointing and ancillary information, using the ‘-fields’ parameter of 
Vanilla. These fields should then be examined closely to find out the max and min of the 
ranges of important parameters such as latitude, longitude and local time. The term 
“ispm[]” should be included to print out the calibrated, interpolated spectrum. 

Questions the user may wish to research at this point include: 
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• Which bodies are visible in the FOV (FOV_TARGETS)? 

• What are the target coordinates of the observation, including latitude 
(LATITUDE_ZPD[5]), longitude (LONGITUDE_ZPD[5]), altitude (ALTITUDE_ZPD[5]) 
etc and are they appropriate? 

• What is the distance to the target (BODY_SPACECRAFT_RANGE), which sets the FOV 
footprint size? 

• What are at the viewing parameters, such as emission angle (EMISSION_ANGLE[5]) 
and phase? (PHASE_ANGLE[5]) 

• Are the Q-points all on the planet, if required (ALL_Q_ON)? 

Finally, the ‘-select’ parameter of Vanilla can be used to narrow the range of data to be 
used in the scientific investigation, using the same fields, with a min and max range 
specified: Vanilla -select “<param> <pmin> <pmax>” (quotes are required). If a specific value 
is required, both min and max ranges have to be set to the same value. 


5.1.3 Bad data check 

There are several methods of checking the data that may be employed: 

• Viewing individual spectra to look for gross distortions/noise. 

• Computing the RMS (standard deviation) of a set of similar spectra: electrical 
interferences (noise spikes) often show up as large departures in the standard 
deviation at certain wavenumbers. See also Section “2.5 Electrical interferences 
(spikes) and other data issues”. 

• Checking that the number of overall raw samples (npts) on each IFGM falls within 3- 
sigma of the mean number (may indicate scan mechanism out of phase lock). 

• Modeling the spectrum with a radiative transfer model. If the model is trusted, then 
it can sometimes show up problems with the data continuum calibration, including 
ripples and sag/droop (unphysical gradients). 


5.1.4 Average spectra 

After the final subset of spectra has been selected, most users will wish to average 
spectra together in one or more ‘bins’ to improve S/N. A common example might be to add 
together all spectra taken during an integration (static dwell observation), for example at a 
particular altitude on Titan’s limb, or for a feature track on Saturn. Flowever, even mapping 
observations are commonly binned, adding together spectra that are recorded at very close 
spatial positions during a slow-moving scan. Often these observations are designed so that 
the movement rate is slow enough that adequate S/N can be obtained post facto by binning 
in spatial bins where the range of FOV center location is equivalent to the FOV size (width). 
For example, a Titan FIRNADMAP is designed so that about 60 5-sec spectra are taken in 300 
seconds, with a FOV movement in this time equal to 4 mrad: the same size as the detector. 
Therefore 60 spectra can be co-added while losing only a factor 2 in spatial resolution. 
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To assess the noise level at each spectral point, the standard deviation (a) may be used 
as a first order-estimate. However, if the spectra being averaged are not homogeneous to a 
high-degree (i.e. if the ‘signal’ is changing from scan to scan as well as the random noise), 
then a will over-estimate the noise level. For this reason it may be desirable to estimate the 
random noise from the a of space calibration spectra, which are taken at very frequent 
intervals (near daily) and included in the database. Note that the expected random noise on 
an average of N spectra is equal to a/V N. 

5.1.5 Field of view modeling 

Before modeling the spectrum, the investigator must also consider whether the field 
sampled by the detector is significantly uniform or non-uniform. By ‘significantly’ we mean 
that the signal non-uniformity is larger than the expected random errors in the spectrum. 
For example, the large FPi footprint (3.9 mrad) when directed to Titan’s limb samples 39 km 
at 10000 km distance (FWHM size), or somewhat less than one atmospheric scale height 
(-50 km). However the time spent at this range is very short, so more typically FPi limb 
observations maybe taken at ranges up to 50000 km, when the FOV is already -200 km in 
size and much larger than a scale height. So in this case, accounting for the non-uniformity 
of the stratified atmosphere within the footprint size when modeling the spectrum would 
be desirable. 

Several journal papers have dealt with this topic in detail, including giving precise 
information on the FOV spatial response functions (Nixon, Infrared Limb Sounding of Titan 
with Cassini CIRS, 2009) (Teanby, 2007), and should be consulted for detailed modeling 
instructions. 


5.1.6 Radiance modeling of the spectrum 

Modeling infrared spectra of planetary atmospheres, surfaces and rings is a large field 
within planetary sciences. Typically this problem can be broken down into two parts (1) the 
forward problem: calculating a reference synthetic spectrum based on a set of model 
parameters, e.g. planetary temperature field and gas composition, and (2) the inverse 
problem, which attempts to adjust the model parameters iteratively, to minimize the model- 
data spectrum mismatch, and thus ‘fit’ the spectrum and derive useful information about 
the emitting object (the model parameters). There are a number of useful textbooks on this 
subject (Goody, 1989) (Hanel, 2003) (Rodgers, 2000). 

In addition, the CIRS data user is advised to carefully read a sample of the already 
published analyses of CIRS data, which give details of several different approaches to 
spectral synthesis. These include the Oxford model NEMESIS, which has been applied to 
model the atmospheres of Jupiter, Saturn and Titan using the correlated-k method Irwin. 
2008}; several line-by-line models used at Observatoire de Paris, Meudon, e.g. AART 
(Coustenis, 2010) ; a second correlated-k model at NASA GSFC (Achterberg, 2011) ; and other 
models. 
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While a full discussion of radiative transfer modeling is beyond the scope of this 

document, a few key pointers specific to modeling the CIRS spectra can be mentioned. 

• The user should be aware of the aliased and unaliased bandpasses of CIRS. See Appendix 
A of fNixon, Titan's Prolific Propane: The Cassini CIRS Perspective. 2Qoq) . 

• The finite field of view of CIRS which convolves the source may need to be modeled 
explicitly, by treating the spectrum as the weighted average of multiple rays emanating 
from different points on the source, if inhomogeneous. See Section “5.1.5 Field of view 
modeling” of this document for further details. 

• The spectral line shape and FWHM should be explicitly modeled. A description of the 
apodization leading to the final shape (Flamming Function) is discussed in Section “3.5 
Post-processing: resampling and apodization”. 


5.2 Worked Example 1: Icy Satellite Science 


5.2.1 Introduction 

This example gives an overview of how the temperature maps of Mimas, used in 
(Howett, 2011) were produced. In this example capitalized words are used to show when 
fields in the pointing, geometry or spectral files are being referred to. FP3 is short for the 
CIRS focal plane 3 detector (similarly for FPi and FP4), NAIF is stands for NASA’s Navigation 
and Ancillary Information Facility and NAIFID is the NAIF’s identification code for solar system 
objects. SPICE is NAIF’s toolkit for interpreting spacecraft data. More details of SPICE can be 
found in (Acton, 1996). RA and DEC are short for Right Ascension and Declination 
respectively. 

5.2.2 Deciding which observations to use 

Icy satellite encounters by Cassini are infrequent enough that the times and dates of the 
encounters are relatively well known. Satellite observations can be searched for using the 
name of the observation since the orbit and target are given in the observation title (e.g. 
CIRS_126MI_X_PRIME will be CIRS PRIME Rev 126 observation of Mimas). SPICE, using the 
Cassini kernels can be used to determine when close encounters between Cassini and a 
target of interest occur. On most encounters some remote sensing occurs, although no CIRS 
data will exist for some gravity or radar focused encounters. At the time of writing Cassini 
SPICE kernels could be downloaded from here: 

http://starbrite.jpl. nasa.gov/pds/viewDataset. jsp?dsid=CO-S/J/E/V-SPICE-6-Vi.o 


Orbit 

Date 

Time 

(UTC) 

Sub-Solar 

Longitude 

Local Time 
Coverage 
between ± 20 
Latitude () 

Range 

(km) 

Average 

Spatial 

Resolution 

(km/pixel) 

126 

13 Feb 2010 

19:10 - 20:30 

0 
r\ 

1 

'Tf- 

kO 

127-264 

37,510 - 65,826 

9.0 


89 


139 

16 Oct 2010 

14:27-16:53 

85-90 

4-130 

73 , 167 - 101,071 

16.0 

144 

31 Jan 2011 

01:30 - 03:05 

34-53 

40 -207 

138,88 - 141,011 

42.2 


Table 20: Details of the CIRS FP3 observations used in this analysis. Local time is defined as the angular rotation of the 
body, in degrees, since local midnight. The names of the observations used are (corresponding from top to bottom in the 
table): CIRS_126MI_FP3 DAYMAPooi_PR1ME, CIRS_139MI_MIMASooi_PRIME and CIRS_144MI_M1MAS001_PRIME. 

In Howett C. e., 2011) data from three flybys of Mimas by Cassini are used, which occurred 
on 13 Feb 2010, 16 Oct 2010 and 31 Jan 2011 and are referred to as Rev. or orbits 126, 139 and 
144 respectively. Details of the observations used, including their full names, are given in the 
paper . However, here we will concentrate on the analysis on Rev 144, since the process was 
the same for all three data sets. The first step in analyzing any data is to download it; in this 
case the pointing, geometry and spectral files (the POI, GEO and ISPM files) were 
downloaded. 

The data was then subseted to include only observations of Mimas by selecting in the 
pointing information those data that had TARGET_BODY for NAIFID of 601 (Mimas), and 
with DET values between 1 and 20 (FP3). Within a single icy satellite encounter there are 
often many types of observations made by CIRS. Usually (where data volume and timing 
allow) an FP3/4 global map is made by slowly scanning across the target at 6 /.i radian per 
second, and an FPi raster scan is made at Nyquist sampling at 40 n radian per second. When 
making a temperature map each of the observations are considered separately, since local 
time changes throughout an encounter smearing surface temperatures. 

During an encounter each of the different observations that occurred can be determined 
by looking at the SCET (CIRS time stamp) versus RA and DEC variations in the pointing 
information. Integrations will have little change in the RA and DEC with SCET, whereas 
global scans with have smoothly varying values. In this example a global FP3 scan is used, 
which is comprised of two limb to limb swaths obtained on the 31 st Jan 2011 between 
02:54:35 and 03:05:30 UTC (i.e. SCET 1296442474 to 1296443130). 

5.2.3 Pointing Adjustments 

Since the icy satellites are much smaller than other targets in the Saturnian system, such 
as Titan or Saturn, pointing errors correlate to greater longitude/latitude errors. Thus, when 
working on icy satellite data pointing tweaks are often required. How these tweaks are 
determined will depend on the target itself, the type of encounter and required output. For 
example on Enceladus it may be acceptable to tweak the pointing so that the warmest 
observed temperatures coincide with the hot sulci on Enceladus’ South Polar Terrain (often 
dubbed tiger stripes). However, since the accurate determination of the boundary of Mimas’ 
thermal anomaly was a key goal of this work, accurate absolute pointing was very 
important. 
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Figure 28: temperatures across Mimas from Rev 144, plotted in 
RA/Dec space with the expected limb location depicted by the 
white circle. 


Currently SPICE is used to check the 
pointing, however other systems exist 
that could perform this function. The 
method used here lines up the predicted 
limb location to that observed in the 
temperature information. This works 
well for a global or limb-to-limb daytime 
observation against dark sky, or for a 
nighttime observation during a Saturn 
transit. It is more limited for a nighttime 
observation against dark sky, where 
there is no large temperature variation 
at the limb or for an observation that 
does not scan off the limb (e.g. a stare). 
However, since FP3 global scans were 
made of Mimas the data technique 
above could be used. 


The best fitting blackbody 
temperature curve was found for each 
spectrum in the designated time frame, and the RA and DEC of each observation read from 
the pointing files. The expected limb location in RA and DEC space at the start of the 
observation was determined using SPICE and the Cassini kernels. Using this information a 
temperature map versus RA/DEC was determined, with the limb prediction overlaid, as 
shown in Fig. 23. The limb position is shown at the start of the scanning time and the fields 
of view have been scaled to this time to account for the changing distance of the 
observation. 


If the fitted blackbody temperatures show a dramatic variation where the limb is 
predicted to occur then no tweaks are required. However, as Figure 23 shows, the predicted 
and observed limbs don’t always overlap. So tweaks in ra and dec space can be performed 
until the expected limb location from SPICE overlays the temperature limb. Using further 
SPICE routines these ra/dec tweaks can be translated into longitude and latitude tweaks. 

5.2.4 Rejection of bad data. 

Occasionally hot pixels make it into the CIRS data set. These individual data points are 
usually extremely obvious in an observation, having much higher temperatures than their 
neighboring, often spatially overlapping, observations. An easy check for these is to plot the 
SCET time against the SPECTRUM radiance (found in the spectral files) at a high signal to 
noise wavenumber (e.g. at -200 cm-i for FPi), or plot this radiance value or its 
corresponding temperature against RA and DEC. As Figure 29 shows, any bad scans are likely 
to be very obvious outliers on these plots and once identified they can be excluded from the 
analysis. 
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5.2.5 Data Analysis 

The first data product produced for icy satellite 
surfaces from CIRS data is often a surface 
temperature map (and is the goal of this example 
analysis). 

In the event that pointing tweaks were not 
required then the data is further limited to ensure 
that the fields of view of the data used lie on the 
surface of the target. This is achieved by selecting 
data that has pointing information with a 
FILLING_FACTOR greater than 0.99, ALL_Q_ON 
equal to 1 and ALTITUDE equal to 0. For FP3 and 
FP4 the field of view boundaries are given in the 
pointing file for each observation by the arrays 
LONGITUDE and LATITUDE. For FPi using the 
information on the FPi field of view in (Flasar, 
2004) or by using SPICE and Cassini SPICE kernels 
the circular field of view can be plotted instead of 
the octagonal one given by the LONGITUDE and 
Figure 29 - Temperature RA/DEC plot for Rev 126 LATITU DE in the pointing files. 

Mimas, showing hot bad pixels in the center of the 

dlsk ' If the pointing has been tweaked then the data 

will still have to be restricted to spectra that have 
their complete field of view on the surface. The latitude and longitude of each observation 
will be determined from the tweak. 


The surface temperature correlated to each of the remaining spectra is determined by 
finding its best fitting black body temperature emission curve. A surface temperature map 
can then be created by combining the observations longitude, latitude and temperature 
information. Where data coverage overlaps the mean temperature in the overlap portion is 
used. For Mimas this process results in the creation of the map shown in Figure 30. 
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Figure 30 - Surface temperature map of Mimas, produced using Rev 144 data. For full description of 
this plot see Figure 2 of Howett et al., 2011. 


5.3 Worked Example 2: Saturn Science 


5.3.1 Overview 

The Cassini Composite Infrared Spectrometer (CIRS) regularly provided 10-1400 cm 1 (7- 
1000 pm) observations of Saturn’s thermal emission throughout the prime and extended 
missions, allowing a detailed characterisation of the seasonal evolution of the temperature 
and composition of a giant planet atmosphere. In this section, we describe the various types 
of Cassini Saturn data and their intended uses, and techniques for analysing spectra from 
the three CIRS focal planes (FPi, FP3 and FP4). In general, FPi is sensitive to tropospheric 
temperatures and para-hydrogen fraction via the collision induced absorption of H 2 and Fie, 
and to the molecular abundances of phosphine (PFH 3 ), ammonia (N F 3 3 ) and methane (CPI 4 ) via 
their rotational transitions. FP3 is generally used for upper tropospheric temperatures in the 
600-680 cm 1 region, and to determine various hydrocarbon species (mainly ethane, C 2 F-I 6 , 
and acetylene, C 2 PI 2 ). Finally, FP4 is used to sound stratospheric temperatures via CPI 4 
emission near 1300 cm 1 , and the tropospheric distribution of PP 1 3 near 1110 cm 1 . The CIRS 
range also encompasses a multitude of trace species, isotopologues and aerosol features. 
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Figure 31: Example of CIRS 2.5 cm-i spectra of Saturn using the mid-infrared focal planes, FP3 and FP4. The upper abscissa 
shows wavelengths in pm. Emission features of CH4 (7.7 pm), ethane (12.3 pm) and acetylene (13.8 pm) are visible, along 
with absorptions due to PH3 and NH3 (8-11 pm) and the collision-induced continuum of H2 (15+ pm). 


5.3.2 Saturn Data Selection 

Nadir observations of Saturn come in four ‘flavors’ (excluding regional maps which are 
sometimes acquired in tandem with other instruments). The type of observations to be 
used is determined by the particular science application (i.e., temperature mapping, 
compositional measurements, etc.), with some examples given below. Particular data types 
can be selected via the RTI (real time interrupt) field, which determines the length of the 
interferograms, and hence the resolution of the spectrum. Data is collected in both the far- 
and mid-infrared focal planes simultaneously, despite the names given to the observations. 

A useful list of all CIRS (indeed all Cassini observations) is contained in the *TOL.ASC files 
contained in the INDEX/ directory of each PDS volume. In this file you can find the start and 
end date/time for each observation, useful when querying the database for a particular 
observational type. These date/times first need to be converted into Cassini SCET time (UT 
seconds since 00:00:00 Jan 1 st 1970). The data can then be extracted using the following 
fields: 

• SCET times (spacecraft elapsed time based on your selected observation) 

• Detector numbers (0 for FPi, 1-20 for FP3, 21-40 for FP4). 

• RTI (specifies the spectral resolution, see below). 

• Target ID = 699, FOV Targets = 64. 

• Select only open-shutter observations (closed-shutter is for calibration only) 

• Select only observations with all q-points on the disc (otherwise one risks 
contamination from deep space). 
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Low Resolution Mapping (FIRMAP) 

These hemispherical mapping observations,_obtained at 15.0 cm’ spectral resolution 
(36/37 RTlX are the only ones that provide complete ‘maps’ of Saturn’s thermal emission. All 
other observations sit and stare at a single latitude as Saturn rotates beneath (useful for 
high resolution spectroscopy, but not for latitudinal mapping). The FIRMAP observations 
require short integration times (because of the low spectral resolution), which allows the 
fields of view to be scanned over Saturn’s disc from a close proximity (17-37 Rs, Saturn radii), 
providing a spatial resolution of 6-8° latitude in the far infrared and 1-2 0 in the mid-infrared. 
The detectors are scanned from north to south along the central meridian as Saturn rotates 
either once or twice on its axis, so a full FIRMAP observation takes 10-20 hours. 

1 

Composition Integration Mid IR Nadir Map Far-IR Nadir Map 

(0.5 cm ') (3 cm 1) (15.5 cm') 

I I I 

45R S 30 Rj 20R S 



Figure 32: diagram showing three types of observations used - COMPSIT 0.5 cm-i ‘sit-and-stare’ integrations at specific 
latitudes, MIRMAP 2.5 cm-i integrations using the same technique, and hemispherical FIRMAPs at 15.0 cm-i resolution 
where the arrays are scanned north-south over the central meridian. 


FIRMAPs are typically used for mapping atmospheric properties such as the tropospheric 
temperature and para-hydrogen distributions, as these affect the shape of the hydrogen- 
helium collisionally-induced continuum between 250 - 680 cm 1 (FPi and FP3), as well as the 
stratospheric temperatures from the methane emission in FP4. They also provide some 
estimate of the distribution of ethane and acetylene from FP3, albeit at a low spectral 
resolution. An example of a Saturn methane emission map from August 2011 is given in 
Figure 33. 
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Figure 33: FIRMAP 15 cm-i observations of Saturn's methane emission associated with the eruption of the storm in 2010 
(Fletcher et al., 2011). Each FIRMAP focuses on a single hemisphere, but captures all longitudes. 


Medium Resolution Composition (MIRMAP) 

MIRMAPS are obtained by pointing at a single latitude and allowing Saturn’s rotation to 
sweep the arrays over all longitudes. The 2.5 cm 1 data (RTI 96/97) provide a useful 
compromise between spatial coverage, spectral resolution and signal to noise 
characteristics, being obtained from an orbital distance of approximately 30RS. Latitudinal 
coverage must be built up slowly over the course of many months of repeated observations 
(see Fletcher et al., 2009, Icarus), and throughout the course of the Cassini prime and 
extended missions the latitudes were covered on several occasions. 

The ability to resolve spectral features of PH 3 , NH 3 , C 2 H 6 and C 2 PI 2 make these observations 
ideal for coaddition and determination of temperatures and abundances at a single latitude 
(assuming zonal homogeneity around a latitude circle). Furthermore, because the FP3 and 
FP4 arrays are oriented north-south during the sit and stare, these observations can be used 
as rudimentary maps of a 10-degree latitudinal swath (e.g., Figure 34, where two MIRMAPs 
are shown next to one another). 
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Figure 34: MIRMAPs can also be used to map out contrasts in thermal emission for a narrow latitudinal width. This figure 
shows two MIRMAPs from July 2011, focusing on methane emission from Saturn’s storm. 
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Medium Resolution Mapping (MIRTMAP) 

From 2009 onwards, a compromise was needed between latitudinal mapping and 
spectral resolution, because the FIRMAPs were becoming less frequent and could not offer 
the possibility of tracking seasonal trends in chemical species. The MIRTMAPs have the 
same spectral resolution as the MIRMAPs (2.5 cm 1 , RTI 96/97), but are stepped across 
latitude bands with the arrays oriented north south. The resulting distribution of 
measurements is shown in Figure 35. 
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Figure 35: Example of MIRTMAP observations at 2.5 cm-i for the northern hemisphere, showing how the arrays are 
stepped from north to south as Saturn rotates beneath the spacecraft. 

The MIRTMAP observations can be used to determine latitudinal distributions of 
temperatures, PH 3 , N H 3 , C 2 H 6 and C 2 H 2 at a single epoch, rather than waiting for a latitudinal 
map to be generated from the sequential MIRMAP observations. 

High Resolution Composition (COMPSIT) 

The highest spectral resolution offered by CIRS are the 0.5 cm 1 observations (RTI 
400/401), acquired in the same fashion as the MIRMAPs (i.e., sitting and staring at a single 
latitude). Some COMPSITs use high emission angles to maximize the stratospheric emission 
from chemical species (i.e., utilizing the limb brightening effect), whereas others use a close- 
to-nadir geometry to maximize tropospheric absorption features. Long integration times 
are required to obtain a sufficient signal-to-noise ratio at this high resolution, so 
interferograms were measured from an orbital distance of 50-60 Rs, leading to an average 
field of view in the far-IR of io° latitude, and 2-4 0 latitude in the mid-IR. Each observation 
focused on a single latitude, so north-south coverage must be acquired over many months. 

COMPSIT observations are typically coadded for a particular latitude circle to enhance 
the signal-to-noise, then used to provide precise determinations of vertical temperature and 
compositional gradients. They are also used to resolve narrow features, such as rotational 
features of methane and HD or emission due to water and trace hydrocarbons. Many 
hundreds of spectra must be coadded to achieve a useable signal to noise, so COMPSITs are 
not suitable for spatial mapping. 

5.3.3 Estimation of Spectral Noise 

Once the required observation is identified, we need to estimate the uncertainty on the 
measurements. This can be done in several steps: 
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[1] Obtain the CIRS single-scan NESR (noise equivalent spectral radiance) and reduce this by 
the square-root of the number of target scans added together (VN). 

[2] Calculate the standard deviation of the individual scans in the average (also useful for 
identifying spectral regions prone to noise spikes). 

[3] Calculate the noise on the deep space scans used for calibration (i.e., the CIRS NESR 
reduced by the number of deep space spectra (VM)). The number of target spectra 
typically exceeds the number of deep space calibration spectra, such that this dominates 
the uncertainty in the radiance. 

The final uncertainty is taken to be the larger of [i]+[3] (added in quadrature) or the 
standard deviation, [2]. Once the spectral uncertainty is defined, one also needs to be aware 
of possible artifacts in the CIRS data that could corrupt interpretations. 

A note on coaddition: In producing coadded spectra, it is wise to use data from a single 
observation type, as shifts in the instrument characteristics can subtly affect the calibration 
techniques used for each spectrum. So, rather than coadding all 401 RTI spectra for the 
month of August 2010, for example, it is better to isolate and individually coadd just the 
dedicated COMPSIT observations during that time period. 

5.3.4 Possible Systematic Errors 

Below is a non-exhaustive list of some of the pitfalls found in Saturn CIRS data during the 
prime and extended mission (see also Section 2.5 Electrical interferences (spikes) and other 
data issues): 

• Clock Spike: FPi 8-Hz spike (the clock spike) at 191 cm 1 in FPi, a result of the RTI 
count (which is i/8 th of a second). This also has low-power harmonics at 16 Hz (382 
cm') and 24 Hz (573 cm '), and in FP3 at 32 (768 cm ') and 40PIZ (960 cm'). 

• BIU Spikes: A repetitive series of 0.5 Hz spikes appear as a comb on the 

interferograms, resulting from the querying of the CIRS Bus Interface Unit (BIU) by 
the spacecraft Command Data Subsystem (CDS) every two seconds (every 16 RTIs) to 
transfer cached CIRS data. The BlU-spikes occur throughout the interferogram with a 
frequency 0.5 Hz, and are Fourier transformed to a regular comb of spikes in 
calibrated spectra with a spacing of approximately 11.9 cm ’. These spikes were 
largely removed from 2006 onwards via new calibration techniques. 

• Sine Wave: A sine ripple occurring on FPi interferograms, suspected to be from a 
thermal source, leads to a large spike in calibrated FPi spectra. Its frequency varies 
between 1 and 18 Hz over time, so its location is not well defined. Plotting the 
variance of the FPi spectra is helpful to avoid regions of the spectrum contaminated 
by the sine wave. 

• Ripples: Discontinuities in the interferograms, particularly at the transition between 
the double-sided region (symmetric about the ZPD) to the unique, single-sided 
region, and at the ends of interferograms, can cause ripples in calibrated spectra. 
These ripples are particularly acute in the low signal regions of FP3 (900-1100 cm'), 
making derivations of the tropospheric gas composition challenging here. 


98 



• Negative FP4 radiances: Thermal variations and laser mode shifts within the 
instrument between the times of deep space observations and the target 
observations can lead to inappropriate calibration references. These manifest 
themselves as baseline shifts for the radiances, and these are particularly visible in 
the 1100-1250 cm 1 region of FP4. If the radiances are largely negative in this region, 
then the spectra are corrupt and a special calibration is required. 

• FPi short-wavenumber issues: The 10-50 cm 1 region of FPi may be prone to 
distortions associated with an as-yet unidentified thermal effect, leading to radiances 
and brightness temperatures which are non-physical. It is recommended to use this 
region in a line-to-continuum mode only. 

5.3.5 Data Analysis 

Once all sources of noise and systematic errors are accounted for, we proceed with 
inversion of the spectra to determine atmospheric temperatures and chemical composition. 
Several inversion techniques are available, including the multi-purpose NEMESIS optimal 
estimation retrieval code developed by Pat Irwin and colleagues at the University of Oxford 
(Irwin, 2008). The first-order product, whether it is from a spatially resolved FIRMAP 
observation, or a co-added MIRMAP/COMPSIT observation at a single latitude, will be the 
temperature field from the troposphere to stratosphere. These are then fixed and used to 
derive molecular abundances of a host of tropospheric and stratospheric species. 

Temperature Retrievals 

Figure 36 shows Saturn’s mean zonal temperature field from the Cassini prime mission, 
2004-2008, assembled by coadding all available FIRMAP observations (see Section 2.3.1). 
Tropospheric temperatures were derived from the 100-600 cm 1 region of FPi, which provide 
a lower spatial resolution dataset that can then be updated with the higher spatial 
resolution of FP3 and FP4. Thus the FP3 600-680 cm' spectrum (collision induced H 2 
absorption, sensitive to 80-200 mbar) and the FP4 1250-1350 cm 1 spectrum (methane 
emission from the 0.5-10 mbar range) are used simultaneously to constrain the upper 
troposphere and stratosphere. Note that a good first-guess stratospheric profile is useful 
when retrieving temperatures from FPi, as the peaks of the collision induced absorption do 
have some sensitivity to lower stratospheric temperatures. 
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Figure 36: Saturn's temperature field derived from a co-addition of all available FIRMAP observations from the prime 
mission. The region between 10-80 mbar is a smooth interpolation between the two regions of CIRS sensitivity. 
Temperatures below 800 mbar are a smooth relaxation back to the a priori profile. See (Fletcher L. N., 2010) for details. 


Compositional Retrievals 

Provided that the temperature profiles are defined for the vertical range of interest, 
gaseous composition can then be adjusted to fit emission and absorption features. Typical 
species of interest include PH 3 (a disequilibrium tracer of tropospheric dynamics), NH 3 (a 
condensable species in the upper troposphere, used to assess cloud condensation 
processes); CH 4 (well-mixed, so can be used to estimate the C/H, D/H and 12 C/ 13 C ratios in 
Saturn’s bulk) and hydrocarbon species (principally C 2 H 6 and C 2 H 2 ). Trace species are also 
present and can be investigated in spectra with many thousands of coadded observations. 
Three different techniques are suggested, depending on the spectral resolution of the 
selected data (low resolutions limit the information content in a spectrum) and the desired 
speed of the spectral inversion: 

[1] Scaling reference profiles: Either obtain a theoretical profile for a species (e.g., the 
saturated vapour profile, or a photochemical profile) and scale it to fit the data. 

[2] Parametric profiles: For some species, define the profile in terms of a deep mole 
fraction, a reference pressure, and a fractional scale height above that level (the ratio 
of the gas scale height to that of the atmosphere). This technique has been 
successfully used for PH 3 on Jupiter and Saturn (Figure 37) (Fletcher L. N., Phosphine 
on Jupiter and Saturn from Cassini/CIRS, 2009) 
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[3] Continuous profiles: For high spectral resolution data where multiple transitions are 
resolved, it is appropriate to derive a continuous vertical profile of species (e.g., 
ethane and acetylene) defined on multiple levels. This technique was used by 
(Howett C. e., 2007). 


Saturn PH 3 Distribution (ppm) 



Planetographic Latitude 

Figure 37: The meridional distribution of PH3 derived from CIRS 1100-1200 cm-i FP4 observations at 2.5 cm-i using a 
parameterised profile (Fletcher L. N., Phosphine on Jupiter and Saturn from Cassini/CIRS, 2009). 
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Further Information Sources 


Hopefully, this document contains the majority of the information you are looking for. If not, 
please contact a CIRS representative listed in “Appendix E. Contact Information”, so long as 
the mission is continuing. After that time, someone at the PDS Atmospheres node may be 
able to help. 
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Appendix 


A. Example PDT Designs For Observations 


In this section we show example screenshots for the design of Cl RS observations using the 
Cassini PDT software. Detailed descriptions are given in Section “2.3 Observation 
Strategies”. 
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Figure 38: COMPSIT with the focal planes on the central meridian. Also shows the deep space position. 
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Figure 39: COMPSIT with FPi offset for high emission angle studies. Also shows deep space position. 



Figure 40: COMPSIT with FP3 offset for high-emission angle studies. Also shows deep space position. 
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Figure 41: MIRMAP on the central meridian. Also shows deep space position. 



Figure 42: MIRTMAP with focal planes at multiple latitude positions on the central meridian. Also shows the deep space 
positions. 
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Figure 43: FIRMAP showing the targeting positions and deep space position. 



Figure 44: FIRMAP showing the focal planes for part of an observation. This shows the slewing across the northern 
hemisphere. 
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Figure 45: LIMBINT showing the entire observation and deep space position. Please note that LIMBMAP's also look like 
this except they have many target latitudes. 
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Figure 46: LIMBINT at altitude position 1. Focal planes are oriented normal to the local latitude limb to sound the first 
altitude position. This will be the same for LIMBMAP's. 
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Figure 47: LIMBINT at the second altitude position. 



Figure 48: NADIROCC shown along with the deep space positions. 
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A. 2 RINGS 
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Figure 49 - A TEMP, CIRS_074RI_TEMPLi3LP00i_PRIME, showing four radial scans taken from a relatively low spacecraft 
emission angle. 
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Figure 50 - 



Figure 51 - A TDIF is composed of two scans: one on the lit side of the rings, the other on the unlit side of the rings but 
with an otherwise comparable geometry. In this case the rings have been scanned at a local hour angle of I2h and i8h 
from both the lit and unlit side of the rings. 
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Figure 52 - Temperatures within Saturn's shadow are measured by scanning the FPi detector parallel to the planet's 
shadow. Each scan maps out locations in the rings that have been within the planet's shadow for the same amount of 
time. 



Figure 53 - Azimuthal scans reveal temperature variations in the course of an orbit around Saturn. Due to RBOT 
considerations, such long duration azimuthal scans are no longer performed. 
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Figure 54 - During POLARIZ observations CIRS integrates on the rings from three distinct orientations in order to look for 
any polarization in the signal originating from the rings. The outer C ring is being targeted in 
CIRS_ii3RI_POLARIZi8i_PRIME. Deep space spectra are also obtained at the same three spacecraft orientations. 
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A.3 TITAN 


Figure 55: T6 FIRLMBT 




Figure 56: T6 FIRLMBAER 



Figure 57: T6 FIRLMBINT 



118 




Figure 58: T6 FIRNADMAP 
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Figure 59: T6 MIRLMBINT 
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Figure 60: T6 MIRLMBMAP 



Figure 61: T6 FIRNADCMP 
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Figure 62: T6 M 1 DIRTMAP 
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A.4 satellites 



Figure 63: Mimas FP3DAYMAP in Rev 126A.4 SATELLITES 
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Figure 64: Enceladus FPiDRKMAP Rev 131 
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Figure 65: Enceladus SECLNX in Rev 8 
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Figure 66: lapetus FP13JVIAPS Rev 00C 
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Figure 67: Hyperion FPiGREEN Rev 15 



Figure 68: Tethys FP1PF3MAP Rev 15 
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Figure 69: Enceladus Collaborative Observation taken in Rev 28. Includes ISS, UVIS, and VIMS pointing as well as CIRS FPi, 
FP3, and FP4. 



Figure 70: Rhea FPISECLNX - Rhea eclipse exit on REV 71 
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Figure 71: Tethys FP13STARE taken in Rev 73 



Figure 72: Enceladus distant FP3/4 map and FPi Stare taken in Rev 74 
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Figure 73: Enceladus FPi Eclipse Scan/Map of the South Pole on Rev 80 
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B. List of major data gaps 


Table 21: list of major data gaps 

Observation name 


CIRS_C22JU_17ATM1X1265_ISS 

Cl RS_000SA_M I RCM PSTA004 _PRI M E 

CIRS_000PH_FP13LTCRV004_PRIME 

CIRS_000RI_URIMPACT001_UVIS 

CIRS_000RI_COMP001_PRIME 

CIRS_00BDI_FP1DRKMAP001_PRIME 

CIRS_006RI_FMONITOR001_PRIME 

CIRS_013HY_GEOLOG009_ISS 

CIRS_013HY_GEOLOG011_ISS 

CIRS_013HY_GEOLOG013_ISS 

CIRS_014TI_M I Dl RTM AP006_PRI ME 

Cl RS_014TI_FIRNADMAP005_ISS 

CIRS_014TI_FIRLMBAER004_UVIS 

Cl RS_0 14TI_M I RLM Bl NT002_PRI M E 

CIRS_014TI_FIRNADCMP002_PRIME 

CIRS_014TI_MIDIRTMAP005_PRIME 

CIRS_0150T_ISSSPACE004_ISS 
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CIRS_023RE_E150PHASE001_VIMS 
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CIRS_023RH_238W159PH001_ISS 

CIRS_023RE_E150PHASE007_VIMS 

CIRS_030IC_NORTH7CAL001_RADAR 

CIRS_030IC_DIORPX001_CDA 
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CIRS_032IC_NORTH4CAL001_RADAR 

CIRS 037IC SOUTH3CAL001 RADAR 
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Time 

Observation End Time 
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11/2/2000 0:05:08 

11/3/2000 3:23:11 

0.00% 

6/8/2004 3:28:00 

6/8/2004 5:28:00 

0.00% 

6/11/2004 4:40:00 

6/11/2004 5:23:00 

0.00% 

7/4/2004 0:31:00 

7/4/2004 2:31:00 

0.00% 

7/4/2004 2:31:00 
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0.00% 

9/7/2005 17:11:57 

9/7/2005 20:11:57 

0.00% 

9/7/2005 20:11:57 

9/8/2005 2:23:49 

0.00% 

9/19/2005 6:25:00 

9/19/2005 6:55:00 

0.00% 

4/30/2006 5:34:15 

4/30/2006 11:58:15 

0.00% 

5/1/2006 23:10:00 

5/2/2006 2:00:00 

0.00% 

5/2/2006 2:40:00 

5/2/2006 3:50:00 

0.00% 

5/2/2006 5:35:00 

5/2/2006 7:35:00 

0.00% 

5/2/2006 9:20:00 
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0.00% 

5/2/2006 13:20:00 

5/2/2006 16:30:00 

0.00% 

10/7/2006 8:00:00 

10/7/2006 9:45:00 

0.00% 

10/12/2006 3:15:00 

10/12/2006 4:45:00 

0.00% 

10/13/2006 23:45:00 

10/14/2006 7:50:00 

0.00% 

10/14/2006 7:50:00 

10/14/2006 9:35:00 

0.00% 

10/28/2006 2:30:00 

10/28/2006 4:30:00 

0.00% 

10/28/2006 4:30:00 

10/28/2006 6:00:00 

0.00% 

11/4/2006 3:12:00 

11/4/2006 5:12:00 

0.00% 

11/4/2006 21:12:00 

11/4/2006 22:42:00 

0.00% 

11/6/2006 20:30:00 

11/6/2006 22:00:00 

0.00% 

11/13/2006 21:14:00 

11/13/2006 22:14:00 

0.00% 

1/15/2007 16:30:00 

1/15/2007 18:00:00 

0.00% 
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CIRS_058RI_RPXMOVIE001_VIMS 
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CIRS 108RF FMOVIEOOl ISS 


1/18/2007 20:20:00 
1/21/2007 1:45:00 
2/21/2007 11:23:25 
2/21/2007 12:12:24 
2/21/2007 14:12:24 
2/21/2007 17:12:24 
2/21/2007 18:12:24 
2/22/2007 14:12:24 
2/22/2007 16:12:24 
2/22/2007 17:12:24 
3/26/2007 1:56:27 
4/21/2007 6:46:00 
4/21/2007 7:50:00 
4/21/2007 13:30:00 
6/8/2007 10:40:00 
9/9/2007 22:05:00 
9/26/2007 21:51:00 
9/27/2007 5:21:00 
9/30/2007 22:30:00 
9/30/2007 23:20:00 
10/2/2007 9:52:43 
10/2/2007 12:42:43 
10/2/2007 13:22:43 
10/6/2007 19:55:00 
11/19/2007 1:08:25 
11/19/2007 1:31:25 
2/10/2008 16:06:00 
4/20/2008 7:32:00 
5/2/2008 4:40:00 
12/17/2008 9:30:00 
3/12/2009 11:45:00 
3/12/2009 20:55:00 
3/13/2009 16:45:00 
3/13/2009 18:45:00 
3/14/2009 1:10:00 
3/15/2009 4:00:00 
3/15/2009 16:20:00 
3/16/2009 4:00:00 
3/16/2009 5:45:00 
3/16/2009 7:45:00 
3/17/2009 3:20:00 
3/17/2009 4:00:00 
4/16/2009 9:13:00 
4/16/2009 10:28:00 


1/18/2007 21:50:00 
1/21/2007 3:15:00 
2/21/2007 12:12:24 
2/21/2007 14:12:24 
2/21/2007 17:12:24 
2/21/2007 18:12:24 
2/21/2007 22:02:24 
2/22/2007 16:12:24 
2/22/2007 17:12:24 
2/23/2007 0:27:24 
3/26/2007 2:38:27 
4/21/2007 7:50:00 
4/21/2007 13:30:00 
4/21/2007 20:45:00 
6/8/2007 17:40:00 
9/9/2007 22:25:00 
9/27/2007 5:21:00 
9/27/2007 9:36:00 
9/30/2007 23:20:00 
10/1/2007 3:20:00 
10/2/2007 12:42:43 
10/2/2007 13:22:43 
10/2/2007 13:42:43 
10/6/2007 20:00:00 
11/19/2007 1:31:25 
11/19/2007 2:47:25 
2/10/2008 19:36:00 
4/20/2008 12:20:00 
5/2/2008 7:40:00 
12/17/2008 11:43:00 
3/12/2009 20:55:00 
3/12/2009 21:55:00 
3/13/2009 18:45:00 
3/14/2009 1:10:00 
3/14/2009 1:50:00 
3/15/2009 16:20:00 
3/15/2009 16:50:00 
3/16/2009 5:15:00 
3/16/2009 7:45:00 
3/16/2009 17:45:00 
3/17/2009 4:00:00 
3/17/2009 17:25:00 
4/16/2009 10:28:00 
4/16/2009 20:33:00 
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4/20/2009 8:50:45 
4/20/2009 9:20:45 
4/20/2009 14:20:45 
4/24/2009 21:00:00 
5/4/2009 0:26:00 
5/4/2009 14:42:00 
8/9/2009 2:01:49 
8/9/2009 4:03:53 
8/9/2009 5:03:53 
8/9/2009 14:29:53 
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8/27/2009 2:45:00 
8/27/2009 6:45:00 
8/27/2009 9:10:00 
8/27/2009 11:30:00 
8/27/2009 15:00:00 
8/27/2009 17:00:00 
8/27/2009 20:45:00 
8/28/2009 0:45:00 
8/28/2009 16:40:00 
8/28/2009 20:40:00 
8/29/2009 2:40:00 
8/29/2009 16:40:00 
8/29/2009 20:25:00 
8/30/2009 16:40:00 
8/31/2009 16:35:00 
8/31/2009 20:20:00 
3/6/2010 19:29:00 
4/18/2010 16:13:00 
7/7/2010 9:22:45 
7/7/2010 14:22:45 
12/4/2010 5:00:00 
12/4/2010 17:00:00 
5/11/2011 9:12:00 
10/2/2007 15:42:43 
7/3/2004 3:15:00 
10/1/2007 18:42:43 
4/19/2009 19:20:45 


4/16/2009 23:32:00 
4/17/2009 20:00:00 
4/17/2009 22:53:00 
4/19/2009 14:20:45 
4/19/2009 15:20:45 
4/19/2009 19:20:45 
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4/20/2009 14:20:45 
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5/4/2009 8:42:00 
5/4/2009 21:36:00 
8/9/2009 4:03:53 
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8/28/2009 6:45:00 
8/28/2009 20:40:00 
8/29/2009 2:40:00 
8/29/2009 7:00:00 
8/29/2009 17:55:00 
8/30/2009 5:49:00 
8/31/2009 6:18:00 
8/31/2009 17:50:00 
9/1/2009 14:15:00 
3/6/2010 20:44:00 
4/18/2010 23:13:00 
7/7/2010 14:22:45 
7/7/2010 17:57:45 
12/4/2010 17:00:00 
12/5/2010 6:00:00 
5/11/2011 10:42:00 
10/2/2007 17:42:43 
7/3/2004 4:00:00 
10/1/2007 19:42:43 
4/19/2009 21:20:45 
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2/1/2009 5:07:00 
4/24/2009 9:55:00 
2/27/2009 4:35:00 
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C. The CIRS FPi Spike Suppression Algorithm: Methodology 


The FPi spike suppression algorithm can be summarized as follows: 

• “Pure” FPi 0.5 Hz and 8.0 Hz spike combs are constructed from “global” large averages 
of deep-space interferograms recorded over a several-month time period. 

• The 0.5 Hz and 8.0 Hz spike combs, the target interferograms (Saturn, Titan, rings, 
satellites, etc.), and the deep-space interferograms are first over-sampled by a factor of 
8. The 8X over-sampled spike combs are then independently adjusted in position and 
amplitude to fit the 0.5 Hz and 8.0 Hz spike patterns in the 8X over-sampled target and 
deep-space interferograms. 

• The 8X over-sampled position- and amplitude-adjusted spike combs are then subtracted 
from the 8X over-sampled target and deep-space interferograms in order to calculate 8X 
over-sampled spike-free interferograms. 

• The 8X over-sampled spike-free target and deep-space interferograms are down- 
sampled to their original number of samples. 

• The down-sampled spike-free target and deep-space interferograms are used to 
calculate spike-free radiometrically-calibrated spectra. 


The logic of the algorithm is illustrated in Figure 74. The “sine wave” suppression is not 
included in the present data base release. It will appear in a future release. 
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From Noise 
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0.5 Hz Spikes Not Suppressed 
8.0 Hz Spikes Not Suppressed 
Sine Wave Not Suppressed 
BIURTI = 66; FIFM ID = 6 



0.5 Hz Spikes Not Suppressed 
8.0 Hz Spikes Not Suppressed 
Sine Wave Not Suppressed 
BIURTI = 77, 88, or 99; FIFM JD = 7 



0.5 Hz Spikes Suppressed 
8.0 Hz Spikes Not Suppressed 
Sine Wave Not Suppressed 
BIURTI = 0-15; FIFMID = 4 



0.5 Hz Spikes Suppressed 
8.0 Hz Spikes Suppressed 
Sine Wave Not Suppressed 
BIURTI = 0-15; FIFM ID = 5 



De-Spiked FIFM 
BIURTI = 0-15 
FIFM ID = 0-3 


Figure 74: CIRS spike suppression algorithm block diagram 
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C.i Calculation of the 0.5 Hz and 8.0 Hz reference waveforms (“combs”) 


Separate FPi 0.5 Hz and 8.0 Hz spike combs are constructed from “global” averages of 
deep-space interferograms. The interferograms in each average are carefully selected to be 
free of external noise effects, such as transient spikes, ZPD shifts, laser mode drifts, and 
velocity variations. The CIRS interferogram noise detector is used in conjunction with 
restricting the lengths of the interferograms to assure optimal quality. The 0.5 Hz and 8.0 
Hz combs are calculated separately from each of 16 average 401 RTI deep-space 
interferograms with the same 0.5 Hz BIURTI offset value. The 0.5 Hz BIURTI offset = the 
number of RTIs that have elapsed from the beginning of sampling until the appearance of 
the first 0.5 Hz spike in the interferogram. An individual interferogram will therefore have 
one of 16 possible BIURTI offsets = 0, 1, 2, 3, ..., 15 RTI. 

The D. C. baseline and the 401 RTI FPi interferogram shape function are first subtracted 
from each of the 16 average deep-space interferograms (BIURTI offsets = 0-15 RTI). The 
shape function is calculated by twice Gaussian-filtering an average of 48,475 401 RTI deep- 
space interferograms that were recorded during the CIRS Saturn mission. Each 
interferogram is then over-sampled by the factor F = 8 and now has 8 times as many points 
as it had originally. The "flat" region of each over-sampled average interferogram is defined 
to begin at the over-sampled sample number N = n*F = 345*8 = 2760 and is far enough 
beyond the interferogram’s ZPD region to be fairly uniform in D. C. level. The 0.5 Hz and 8.0 
Hz spike patterns in the “flat” region are now separated to form the respective spike combs. 
Each comb is then extrapolated from the “flat” region to cover the entire interferogram 
range, plus ~ 410 samples (after down-sampling) beyond the normal length of a 401 RTI 
interferogram. The extra length ensures that a 0.5 Hz comb of any BIURTI offset may be 
fitted to an interferogram of the same or any other BIURTI offset for the spike suppression. 
The fifteen 0.5 Hz combs with BIURTI offsets = 1, 2, 3, ..., 15 RTI are next shifted to coincide 
with the comb with BIURTI offset = 0, and the sixteen coincident combs are averaged. The 
sixteen 8.0 Hz combs are averaged without being shifted. Next, each 0.5 Hz spike in the 
averaged extrapolated comb is replaced by the average of the 0.5 Hz spikes in the “flat” 
region. The “global” 0.5 Hz and 8.0 Hz spike combs that were used to de-spike all of the 
spectra in the current PDS database were calculated from the following deep-space 
interferograms: 401 RTI, interferogram length = 2770 - 2780 samples, DSCAL = 1, BIURTI 
offset = 0-15 RTI, 8,758 interferograms averaged from 10/01/2005 - 02/28/2006, CIRS 
reference laser mode = 02. 

C.2 Calculation of de-spiked filtered interferograms (FlFMs) 


The D. C. baseline and the FPi interferogram shape function are first subtracted from a 
single or average interferogram that is to be de-spiked. If an average interferogram is to be 
de-spiked, then each interferogram in the average must have the same 0.5 Hz BIURTI offset 
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value. The interferogram is then over-sampled by the factor F = 8. The 0.5 Hz and 8.0 Hz 
Spike combs are already over-sampled. 

The BIURTI offset value of the over-sampled interferogram that is to be de-spiked is first 
determined by a cross-correlation between the interferogram’s 0.5 Hz spike pattern and the 
spike pattern in the over-sampled 0.5 Hz comb with BIURTI offset = 0. The 0.5 Hz comb is 
then shifted by the number of samples required to make its 0.5 Hz spikes coincide with the 
spikes in the interferogram. A similar cross correlation and shift of the 8.0 Hz spike pattern 
is subsequently performed. 

"Flat" regions of the shifted over-sampled 0.5 Hz and 8.0 Hz spike combs and the 
interferogram are now defined. The "flat" regions begin at the over-sampled sample 
number N = n*F = 345*8 = 2760 and are far enough beyond the interferogram’s ZPD region 
to be fairly uniform in D. C. level. 

The “flat” regions of the shifted 0.5 Hz and 8.0 Hz over-sampled spike combs are 
independently adjusted in position and amplitude to fit the 0.5 Hz and 8.0 Hz spike patterns 
in the “flat” regions of the over-sampled interferogram. In the course of this process, which 
is performed iteratively, the 0.5 Hz and 8.0 Hz spike patterns in the interferogram are 
separated from each other in order to be independently fitted to the 0.5 Hz and 8.0 Hz spike 
combs, respectively. It was found that only two iterations are required to closely match the 
positions and amplitudes of the 0.5 Hz and 8.0 Hz spikes in the combs to the corresponding 
spikes in the interferogram. If a single interferogram is being de-spiked, some of the 0.5 Hz 
spikes may be missing or greatly reduced in amplitude due to noise. Therefore, the 
correlation is verified by calculating the sample numbers and negative amplitudes of the 
interferogram’s 0.5 Hz spikes that occur within ± F/2 = ± 8/2 = ± 4 samples of the position of 
each 0.5 Hz spike in the shifted 0.5 Hz comb. A minimum of twelve 0.5 Hz spikes with 
amplitudes < 4.0 standard deviations from the mean must be present in the “flat” region of 
a 400 or 401 RTI interferogram in order for the correlation to be considered successful. The 
adjusted amplitude of the 0.5 Hz comb is calculated from the ratio of the average spike 
amplitudes in the shifted comb and the interferogram. The adjusted amplitude of the 8.0 Hz 
comb is calculated from the square root of the ratio of the 8.0 Hz power spectra of the 
comb and interferogram, including the 16.0 Hz and 24.0 Hz harmonics. 

The “flat” regions of the over-sampled, shifted, and amplitude-adjusted 0.5 Hz and 8.0 
Hz spike combs are extrapolated to the samples range of the entire over-sampled 
interferogram = ZPD region + “flat” region. The entire over-sampled, shifted, and amplitude- 
adjusted 0.5 Hz and 8.0 Hz spike combs (ZPD region + “flat” region) are then subtracted 
from the entire over-sampled interferogram. Finally, the de-spiked over-sampled filtered 
interferogram (FIFM) is down-sampled to its original length and is now ready to be included 
in the calibration process. 
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D. Effective Integration Times 


The method for calculating an effective integration time for a scan of a given spectral 
resolution was developed by L. Herath. CIRS interferograms are not purely one-sided or two- 
sided, since data sampling is initiated ~ 1.5 seconds before the carriage reaches the ZPD 
point (x = L,). The data sampled before ZPD will be averaged with the corresponding data 
sampled after ZPD in the ~ symmetrical two-sided region centered about ZPD. The total 
number of data samples acquired after ZPD determines the spectral resolution. The pre- and 
post-ZPD interferogram regions are illustrated in Figure 75, which is an average of FPi deep 
space interferograms recorded December 6 - 9, 2000. 



All interferogram data do not contribute equally to the spectral noise when apodization 
is applied. Data near ZPD (x = L,) contribute the most to the noise, whereas data at 
maximum path difference (x = L, + L 2 ) contribute very little. Conversely, adding integration 
time near ZPD contributes more to reducing the spectral noise than does adding integration 
time near maximum path difference. Therefore, the extra 1.5 seconds of integration time 
before the carriage reaches the ZPD point contributes more than the interferogram average 
to reducing the spectral noise and results in an effective integration time = tsff that is slightly 
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greater than the total scan time = t, + t 2 . This effective integration time is used to calculate 
NESRs for single and averaged scans. The NESRs thus calculated have been verified by 
comparisons with standard deviations of large averages of deep-space spectra. 


The effect of apodization is expressed as a factor q, with the spectral noise proportional to 

Vq: 

dx (°- 1 ) 


f(x) = apodization function 
L, = scan length before ZPD (cm) 

L 2 = scan length from ZPD to the end of scan (cm) 
Equation (D.i) may be expressed in terms of time as: 


<7 = 


1 

*\ +t 2 



A(t) = apodization function 


-( 0 . 2 ) 


t, = scan time before ZPD (sec) 

t 2 = scan time from ZPD to the end of scan (sec) 


To evaluate the effective integration time, we shall compare q calculated by Equation (D.3) 
to that of a one-sided scan = q 0 , beginning at the ZPD point and ending at t 2 : 

q 0 = -)\mU t (D-4) 

‘2 0 

We now average the pre-ZPD side of the interferogram from 0 to t, with its counterpart on 
the post-ZPD side from t, to 2t„ thereby reducing the noise in the interval 0 to t, by V 2 . q 
then becomes: 


1 



/|(i/V2)/(o| 2 *+/|/w| 2 * 


The effective integration time is then: 


t E fr =t 2(q 0 /qi) 


_(D-5) 


(D.6) 
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For the Cl RS instrument a Hamming apodization function is used to calculate t E ff: 


m 


= 0.54 + 0.46* cos 



tMax = maximum carriage travel time from ZPD to the end of scan = t 2 (sec) 


(D-7) 


It was desired to derive an empirical relationship between t E ff and apodized spectral 
resolution = R Apod . To this end, tEff was calculated individually for 14 different resolutions 
from 0.50 cm 1 to 20.291 cm 1 using Equations (D.6) and (D.7). The FPi pre-ZPD scan time t, = 
1.506 seconds was used in the calculations. The 14 pairs of {t E ff, R Apod } values were then 
fitted to an inverse 5 th order polynomial in order to empirically determine tEff as a function of 

^ Apod : 




a\ a2 a 3 

= C/0 H + y + y + 

D D z D 3 

^ Apod ^ ' Apod ^ Apod 


a4 


v Apod 


a5 


v Apod 


(11) 


The inverse polynomial fits are illustrated in Figure 76, and the coefficients are listed in Table 
22. Values of tEff for selected spectral resolutions are included in Table 17. 
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t Eff = Effective Integration Time 
Inverse 5th Order Polynomial Fit 


Inverse 5th Order Polynomial Fit Parameters: 

a0 = 0.664011 

al =41.177183 

a2 = -51.049473 

a3 = 60.293167 

a4 = -30.488605 

a5 = 5.515696 

R = Apodized Resolution = 1/2L 2 

t eff = aO + (al/R) + (a2/R**2) + (a3/R**3) 

+ (a4/R**4) + (a5/R**5) 



0 1 2 3 4 5 6 7 8 9 10 1 1 12 13 14 15 16 17 18 19 20 21 22 

Apodized Resolution (cm' 1 ) 


Figure 76: inverse polynomial fit of scan time to spectral resolution 

| Table 22 : Coefficients for the inverse 5th order polynomial fit to effective integration time as a function of spectral 
resolution 


Coefficient 

Value 

ao 

0.664011 

ai 

41.177183 

a2 

-51.049473 

33 

60.293167 

34 

-30.488605 

35 

5.515696 
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Nicolas Gorius 

NASA Goddard Space Flight Center Code 693, Greenbelt, MD 20771, USA 
301-286-0021; Nicolas.gorius@nasa.gov 

PDS Atmospheres Node (master PDS database copy) 

Lyle R. Fluber 

Astronomy Department, New Mexico State University, P. O. Box 30001, Dept 4500. 
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